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A growing  class  of  rare-earth  and  actinide  compounds  containing  localized 
magnetic  moments  exhibit  anomalous  “heavy  non-Fermi-liquid  ’ behavior  charac- 
terized by  a T In  T low-temperature  variation  of  the  heat  capacity  and  a resistivity 
which  approaches  its  zero-temperature  limit  as  Tu  where  1 < v < 2.  We  propose  a 
mechanism  for  such  non-Fermi-liquid  behavior  based  on  “magnetic  frustration,”  a 
feature  of  many  rare-earth  and  actinide  systems.  Using  numerical  renormalization- 
group  methods,  we  study  a Rondo  model  of  three  frustrated  spin-  \ impurities  cou- 
pled to  a nonmagnetic  conduction  band.  Two  distinct  non- Fermi-liquid  regimes 
are  found,  the  first  related  through  an  interchange  of  spin  and  charge  degrees  of 
freedom  to  the  low-temperature  limit  of  the  two-channel  Rondo  model  for  a single 
spin- 1 impurity,  while  the  second  regime  appears  to  be  completely  novel. 
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CHAPTER  1 
INTRODUCTION 


Magnetic  impurities  have  fascinated  physicists  since  the  early  1930s  when  ex- 
perimentalists noticed  a low-temperature  minimum  in  the  electrical  resistivity  of 
noble  metals  which  contained  very  low  concentrations  of  transition  metal  ions  [26] . 
These  materials  were  also  observed  to  exhibit  a Curie  law  1/T  or  Curie- Weiss  law 
l/(T  + 9)  contribution  to  the  susceptibility  at  high  temperatures  T,  indicating  the 
transition  metal  ions  form  free  or  interacting  localized  magnetic  degrees  of  free- 
dom (local  moments).  The  first  attempts  to  explain  the  minimum  in  resistivity 
described  the  impurities  as  centers  of  potential  scattering,  but  were  able  only  to 
generate  a resistivity  which  increased  monotonically  with  temperature.  By  the 
1950s  it  became  clear  the  magnetic  properties  of  the  impurities  in  these  systems 
needed  to  be  taken  into  account.  This  was  first  done  successfully  by  J.  Kondo  in 
1964  [34].  Kondo  showed  the  interaction  of  a local  moment  with  the  spin  degrees 
of  freedom  of  a nonmagnetic  conduction  band  indeed  can  cause  the  resistivity  to 
rise  with  decreasing  T.  His  work  confirmed  the  importance  of  the  magnetic  prop- 
erties of  the  impurities.  It  also  demonstrated  how  interactions  with  a dynamical 
impurity  can  produce  correlations  among  the  conduction  electrons  of  an  otherwise 
noninteracting  conduction  band. 

In  the  late  1970s  and  early  1980s  great  attention  focused  on  a class  of  “heavy 
fermion”  compounds  containing  rare-earth  or  actinide  ions  with  local  moments 
formed  by  4/  or  5/  electrons,  respectively. 
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These  systems  can  be  described  within  a Fermi-liquid  picture,  where  the  low- 
energy  excitations  behave  like  noninteracting  electrons,  except  that  various  prop- 
erties of  the  electrons  are  renormalized  by  strong  correlations.  The  best  example 
of  how  the  low-temperature  properties  of  these  heavy  fermion  compounds  differ 
from  those  of  a conventional  metal  lies  in  the  heat  capacity.  In  a noninteracting 
metal  the  electronic  heat  capacity  takes  the  form, 

c(T)  = ir , (i-i) 

where  7 is  a temperature-independent  constant.  In  the  heavy  fermion  compounds, 
the  low-temperature  linear  coefficient  of  the  specific  heat,  7 = C/T,  tends  to  be 
hundreds  of  times  that  of  normal  metals.  This  implies  a similar  enhancement  of 
the  effective  mass  of  the  conduction  electrons  in  these  compounds.  In  general 
one  categorizes  a compound  with  a linear  specific-heat  coefficient  greater  than  400 
(mJ/moleK2)  as  a heavy  fermion  compound,  though  rare-earth  compounds  exhibit 
values  as  high  as  1620  in  CeAl3,  1600  in  CeCu6,  and  1100  in  CeCu2Si2,  while  the 
values  observed  in  some  actinides  include  1100  in  UBei3,  420  in  UPt3,  and  400  in 
U2Zn17  [24], 

Due  to  the  high  concentration  of  local  moments  in  heavy  fermion  compounds, 
interactions  between  moments  can  play  a role.  Some  materials,  such  as  U2Zni7  and 
UCdn,  exhibit  antiferromagnetic  ordering,  while  others,  such  as  CeCu6  and  possi- 
bly CeAl3,  remain  paramagnetic  down  to  the  lowest  temperatures  investigated  [24], 
A few  materials,  such  as  Ce2Ru2Ge2  [27],  exhibit  ferromagnetic  ordering.  More 
surprisingly,  investigators  have  observed  superconductivity  below  one  kelvin  in 
CeCu2Si2  [60],  UBei3  [51],  UPt3  [61],  and  URu2Si2  [52],  The  appearance  of  super- 
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conductivity  in  heavy  fermion  systems  is  extremely  anomalous,  since  the  presence 
of  local  moments  is  normally  thought  to  lead  to  pair-breaking  and  the  destruction 
of  superconductivity. 

Since  the  mid-1980s  much  experimental  effort  has  gone  into  the  identifica- 
tion and  classification  of  a wide  variety  of  systems  which  show  departures  from  the 
Fermi-liquid  picture.  Examples  include  the  normal  state  of  some  high-Tc  super- 
conductors [42],  one-dimensional  conductors  [70],  edge  states  in  the  quantum  Hall 
effect  [70],  low-dimensional  semiconductor  structures  (“quantum  dots”)  [19],  and 
tunneling  through  metallic  nanoconstrictions  [21].  An  increasing  number  of  rare- 
earth  and  actinide  compounds  have  also  been  shown  to  exhibit  non- Fermi-liquid 
behavior.  Generally,  these  systems  exhibit  a heat  capacity  C oc  — TlnT,  while  the 
resistivity  takes  the  form 


p(T)=po  + AT",  (1.2) 

with  l < v <2  (whereas  u - 2 in  a Fermi  liquid).  A variety  of  mechanisms  have 
been  proposed  for  non- Fermi- liquid  behavior  in  the  heavy  fermion  compounds:  the 
two-channel  Hondo  effect  [19],  the  presence  of  disorder  [44],  and  proximity  to  a 
quantum  critical  point  [62]. 

The  single-impurity  two-channel  Hondo  model  can  reproduce  many  prop- 
erties of  a number  of  U-  and  Ce-based  systems.  However  there  are  two  strong 
objections  to  this  model:  fine-tuning  of  model  parameters  is  required  to  ensure 
non- Fermi-liquid  behavior,  and  interactions  between  impurities  destroy  such  be- 
havior in  a two- impurity  version  of  this  model  [29]. 
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This  thesis  investigates  an  alternative  route  to  non- Fermi-liquid  behavior 
which  retains  the  framework  of  magnetic-impurity  physics,  but  which  avoids  pa- 
rameter tuning  and  which  incorporates  an  important  feature  of  the  magnetic  in- 
teractions in  many  heavy  non-Fermi-liquid  systems,  namely,  magnetic  frustration. 
Our  study  of  a three-impurity  Kondo  Hamiltonian,  the  simplest  possible  model 
which  incorporates  frustration,  reveals  a novel  form  of  non- Fermi-liquid  behavior 
which  is  much  more  robust  than  that  obtained  from  the  two-channel  Kondo  model 
and  related  models. 

We  note  that  similar  considerations  have  been  applied  in  a effort  to  explain 
the  non- Fermi-liquid  properties  of  metallic  nanoconstrictions  [19].  It  has  been  pro- 
posed that  such  systems  can  be  described  in  terms  of  electron-assisted  tunneling  of 
local  defects  [53].  The  simplest  model,  involving  tunneling  backwards  and  forwards 
between  two  defect  levels,  can  be  mapped  directly  onto  the  single-impurity,  two- 
channel  Kondo  Hamiltonian  [19].  Once  again,  fine-tuning  of  model  parameters  is 
required  to  reproduce  the  experimental  properties.  In  work  carried  out  in  parallel 
with  (and  independently  of)  our  study,  Moustakis  and  Fisher  [45]  have  shown  that 
two- channel-  Kondo  properties  can  be  reproduced  without  parameter  tuning  in  a 
model  of  tunneling  between  three  defect  levels.  It  should  be  emphasized  that  Mous- 
takis and  Fisher’s  model  is  considerably  simpler  than  the  three-impurity  Kondo 
Hamiltonian  because  the  former  has  no  impurity  degrees  of  freedom.  Nonetheless, 
the  similar  role  played  by  geometry  in  the  two  cases  is  intriguing,  and  may  point 
to  a more  general  pattern. 

In  the  remainder  of  this  chapter,  we  attempt  to  develop  an  acquaintance  in 
the  reader  with  the  issues  relevant  to  our  work.  We  first  review  the  motivation  for 
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studying  the  effect  of  magnetic  impurities  on  a nonmagnetic  conduction  band.  We 
introduce  the  Hamiltonian  Kondo  applied  to  such  a system,  discuss  the  techniques 
used  to  solve  this  problem,  and  describe  the  low-temperature  physics  contained 
in  its  solution.  We  then  motivate  multiple-impurity  Kondo  models.  We  focus 
on  the  two-impurity  model  and  give  a detailed  description  of  its  low-temperature 
physics.  We  then  turn  to  a discussion  of  non- Fermi-liquid  behavior  in  rare-earth 
and  actinide  systems  and  present  the  models  which  thus  far  have  been  associated 
with  these  systems.  We  review  the  concept  of  magnetic  frustration  and  explain 
why  we  believe  it  could  produce  non-Fermi-liquid  behavior  in  materials  which 
incorporate  local  moments.  We  introduce  the  three-impurity  Kondo  model  we 
have  studied  and  discuss  the  behaviors  one  might  expect  from  this  model.  We 
conclude  with  an  outline  of  the  rest  of  the  dissertation. 

1.1  The  Resistance  Minimum 

As  mentioned  above,  a minimum  in  the  temperature  dependence  of  the  elec- 
trical resistivity  of  not-very-pure  silver  and  gold  [26]  puzzled  physicists  for  thirty 
years.  The  presence  of  a Curie- Weiss  law  contribution  to  the  susceptibility,  the 
linear  dependence  on  impurity  concentration  of  the  depth  of  the  minimum,  and 
observations  that  the  ratio  of  the  T = 0 resistivity  to  the  minimum  resistivity  was 
independent  of  impurity  concentration  [57]  led  J.  Kondo  to  conclude  [34]  interac- 
tions involving  the  magnetic  degrees  of  freedom  of  individual  impurities  and  the 
conduction  band  are  responsible  for  the  resistivity  minimum.  Thus  to  reproduce 
the  minimum  in  resistivity  he  used  the  s-d,  model,  which  describes  a Heisenberg 
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exchange  interaction  between  a noninteracting  conduction  band  and  a local  mo- 
ment. 

The  s-d  model  was  originally  formulated  by  C.  Zener  to  describe  ferromag- 
netism in  nonmagnetic  metals  due  to  the  presence  of  transition  metal  ions  [76].  It 
describes  an  exchange  interaction  between  the  local  moments  of  transition  metal 
ions  and  the  metallic  host  by  analogy  to  magnetic  insulators.  In  the  s-d  model  the 
coupling  between  the  s electrons  of  the  host  and  the  d electrons  of  the  transition 
metal  ion  takes  the  form 


V3d  = -JS-s  = -J  Szsz  + US+s~  + S~s+) 


(1.3) 


Here  S±  = Sx  ± iSy  and  the  interaction  energy  depends  on  the  relative  orientation 
of  the  impurity  spin  S and  the  spin  s of  the  conduction  electrons  at  the  impurity 
site.  The  first  term  on  the  right  hand  side  of  Eq.  (1.3),  which  depends  only  on  the 
2-components  of  the  spins,  describes  a static  interaction.  However,  the  second  and 
third  terms  represent  dynamical  ( “spin- flip” ) interactions  which  can  change  the 
values  of  Sz  and  sz  while  conserving  their  combined  spin.  In  all  cases  of  interest 
the  coupling  J is  antiferromagnetic  (negative). 

The  s-d  model  Kondo  used  to  describe  the  effect  of  dilute  magnetic  impu- 
rities on  nonmagnetic  metals  incorporated  a single  spin-^  impurity  coupled  to  a 
noninteracting  conduction  band,  isotropic  in  momentum  space,  with  an  energy 
range  ±D  about  the  Fermi  energy, 


Hsd  — • S. 

k.tr.cr' 


(1.4) 
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The  conduction-electron  operators  ck(7  annihilate  Bloch  states  of  energy  ek  with 
momentum  k and  spin  2-component  a.  The  conduction-electron  operators  satisfy 
fermionic  anticommutation  relations, 


where  {A,  B)  = AB  + BA.  One  forms  the  conduction-electron  spin  operators  from 
Pauli  matrices  a, 


Kondo  calculated  the  impurity  contribution  to  the  resistivity  as  a function 
of  temperature  using  a third-order  perturbation  theory  expansion  in  the  product 
of  the  density  of  states  at  the  Fermi  energy  q(cf)  and  the  exchange  coupling 
parameter  J.  He  then  added  the  resistivity  due  to  phonons  in  the  host  metal  to 
find 


The  first  term  in  this  expansion  represents  the  phonon  contribution.  The  second, 
temperature-independent  term  describes  J-dependent  potential  scattering  off  the 
impurity,  with  cimp  the  impurity  concentration  and 


{cL  ckV}  = $(k  - k')<W, 


(1.5) 


(1.6) 


p(T)  = aT 5 + cimpb[  1 - AJQo{eF)  In (kBT/2D)\. 


(1.7) 


97r  mJ2 
8 e2heF 


(1.8) 
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where  m and  e are  the  mass  and  charge  of  an  electron  and  h is  Planck’s  constant. 
The  third  term  represents  the  contribution  of  multiple  “spin-flip”  interactions  as 
described  in  Eq.  (1.3). 

Equation  (1.7)  predicts  a minimum  in  the  resistivity  at  a temperature 


T = 

min 


Ag(eF)Jb 
5 a 


^ imp 


1/5 


oc  c. 


1/5 
imp ’ 


(1.9) 


as  observed  experimentally.  Thus  Kondo  demonstrated  the  presence  of  magnetic 
impurities  could  account  for  the  observed  low-temperature  behavior  in  these  sys- 


tems. 


1.2  The  Kondo  Problem 

As  the  temperature  decreases,  the  third-order  term  in  Eq.  (1.7)  diverges. 
Perturbation  theory  breaks  down  when  the  third-order  term  equals  the  second- 
order  term,  at  a temperature  termed  the  Kondo  temperature  Tk'. 

kBTK  ~ Dexp[l/Q(eF)J],  (110) 

where  kB  is  Boltzmann’s  constant.  The  divergences  at  TF  reflect  the  presence  of 
many-body  processes  which  involve  the  entire  conduction-electron  sea. 

The  breakdown  of  perturbation  theory  in  Kondo’s  model  became  known  as 
the  Kondo  problem.  Early  attempts  to  overcome  this  problem  focused  on  removing 
the  logarithmic  divergence  in  perturbation  theory  due  to  the  integral 

[D  dE=  ( D \ 

JksT  E n\kBTj  ' 


(1.11) 
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These  included  work  by  Abrikosov  [1]  to  incorporate  the  effect  of  higher-order  terms 
through  infinite  sums,  Suhl  [63,  65,  67,  66]  on  an  application  of  the  Chew- Low  the- 
ory of  meson-nucleon  scattering  to  the  Kondo  problem,  and  Nagaoka,  Hamann, 
and  Bloomfield  [47,  28,  18]  on  equations  derived  from  double-time  Green’s  func- 
tions. For  the  work  we  present  in  this  dissertation,  a series  of  papers  by  Anderson 
and  others  [8, 10,  11,  75]  — which  tried  to  solve  the  Kondo  problem  by  mapping  the 
s-d  model  to  a one-dimensional  Coulomb  gas  — are  most  important.  Anderson’s 
subsequent,  more  direct  approach  to  these  calculations  [9]  provides  the  foundation 
for  the  renormalization-group  technique  we  employ  to  solve  our  Kondo  model. 

Anderson’s  “poor  man’s  scaling  approach”  [9]  attempts  to  overcome  the  loga- 
rithmic divergence  in  properties  of  the  s-d  model  as  the  temperature  is  decreased  by 
incorporating  the  contributions  of  processes  involving  intermediate  states  near  the 
band  edge  into  the  value  of  the  exchange  coupling  parameter.  Anderson  calculated 
a second-order  scaling  equation  for  the  renormalized  exchange  by  requiring  that 
the  resolvent  Green’s  function  remains  invariant  as  the  bandwidth  is  reduced.  This 
scaling  equation  indicates  the  antiferromagnetic  exchange  interaction  increases  in 
strength  as  one  decreases  the  effective  bandwidth: 


dJ 

d\nD 


ic'i'O')-/2. 


(1.12) 


Since  decreasing  temperature  decreases  the  range  of  allowed  excitations  (i.e.,  the 
effective  bandwidth)  through  the  Boltzmann  distribution,  Eq.  (1.12)  indicates  the 
magnitude  of  the  antiferromagnetic  exchange  coupling  increases  as  the  temperature 
decreases. 
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Anderson’s  perturbative  technique  cannot  describe  the  s-d  model  at  tempera- 
tures below  the  Kondo  temperature,  since  the  strength  of  the  exchange  interaction 
becomes  large  at  TK ■ However  it  was  argued  by  analogy  to  other  systems  that  this 
scaling  equation  should  hold  down  to  T = 0 and  the  exchange  coupling  parameter 
should  go  to  strong  coupling  ( J — ► -oo).  The  low-temperature  behavior  (T  < TK) 
of  the  Kondo  model  was,  therefore,  understood  to  be  due  to  the  formation  of  a spin- 
singlet  from  the  impurity  and  conduction  electrons  localized  about  the  impurity. 
This  agreed  with  observations  of  a constant  Pauli  paramagnetic  contribution  to 
the  susceptibility  for  T — > 0 [69],  and  a power-law  temperature  dependence  of  the 
resistivity  below  TK,  indicating  that  the  magnetic  impurities  become  structureless 
potential  scatterers  at  T — 0. 

In  1974  and  1975  K.  G.  Wilson  published  quantitative  results  for  the  low- 
temperature  behavior  of  the  Kondo  model  [73,  72],  His  renormalization-group 
technique  provided  a nonperturbative  solution  of  the  Kondo  Hamiltonian  which 
avoided  the  logarithmic  divergence  entirely.  Wilson  found  the  low-temperature 
(T  < Tk)  behavior  of  the  Kondo  model  resembles  that  of  a noninteracting  con- 
duction band  subject  to  a phase  shift. 

The  phase  shift  in  the  spectra  Wilson  generated  is  analogous  to  that  cal- 
culated for  electronic  wave  functions  which  undergo  potential  scattering.  One 
associates  the  phase  shift  of  electrons  which  scatter  off  a center  of  potential  scat- 
tering with  the  number  of  electrons  the  potential  localizes  around  the  center.  The 
Friedel  sum  rule  codifies  the  relation  between  the  number  of  localized  conduction 
electrons  and  the  phase  shift  at  the  Fermi  surface  [22], 
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riimp  = 2 [ F deAg(e)  = 2rj(eF)/ir.  (1.13) 

Here  the  number  of  electrons  nirnp  localized  around  an  impurity  site  is  related  to 
the  change  in  the  density  of  states  A g due  to  the  presence  of  an  impurity.  This  in 
turn  relates  to  the  phase  shift  rj  at  the  Fermi  surface. 

The  spectra  Wilson  calculated  for  the  Kondo  model  indicated  a 7r/2  phase 
shift  occurs  in  conduction  electrons  which  scatter  off  the  impurity,  i.e. , nimp  — 1. 
This  and  the  fact  the  ground  state  of  the  model  has  spin  zero  confirmed  conduction 
electrons  quench  the  dynamical  (spin)  degrees  of  freedom  of  the  impurity  at  T = 0 
through  the  localization  of  a single  conduction  electron,  which  forms  a spin-singlet 
virtual  bound  state  with  the  impurity.  The  screening  of  an  impurity’s  spin  through 
the  localization  of  conduction  electrons  around  the  impurity  site  is  now  known  as 
the  Kondo  effect. 

From  the  spectra  he  generated,  Wilson  calculated  the  T = 0 impurity  con- 
tribution to  the  susceptibility  and  the  specific  heat.  Since  the  impurity  spin  is 
quenched,  it  does  not  contribute  directly  to  the  T = 0 susceptibility.  However,  the 
correlations  induced  in  the  conduction  band  renormalize  the  Pauli  paramagnetism 
of  the  band  electrons,  producing  an  additional  susceptibility 

XimP(0)  = (ggB)2w/(ikBTK).  (1.14) 

Here  is  the  Bohr  magneton,  g is  the  (/-factor  for  the  impurity,  and  w is  Wil- 
son’s number,  w = ec+1/4/7r3/2  with  C — 0.5772.  Wilson  found  the  impurity 
contribution  to  the  linear  specific-heat  coefficient  is 
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7 imp  = ^wIcb/STk- 


(1.15) 


From  these  calculations  he  derived  the  value  of  the  Sommerfeld  or  Wilson  ratio 


Rw  = 


4?r  kg  Ximp 


(1.16) 


3 'limp 

which  equals  one  for  a noninteracting  system.  Wilson  found  the  presence  of  cor- 
relations among  the  conduction  electrons  renormalizes  the  Sommerfeld  ratio  so  it 
equals  2 in  the  Kondo  model. 

The  center  of  potential  scattering  formed  by  the  quenched  impurity  has  a res- 
onance at  the  Fermi  surface  with  a width  proportional  to  the  Kondo  temperature. 
Therefore,  the  calculated  impurity  contribution  to  the  low-temperature  resistivity 
goes  as  T2, 


Pimp(T)  - b 


7 r4rc2 
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(1.17) 


We  define  b in  Eq.  (1.8). 

Above  and  throughout  the  rest  of  this  work  we  focus  on  the  impurity  contri- 
bution to  each  property  of  a system,  defined  to  be  the  total  value  of  the  property 
in  the  presence  of  the  impurity  minus  the  total  value  of  the  property  in  the  absence 
of  the  impurity. 


13 


1.3  Interactions  Among  Magnetic  Impurities 

As  we  outlined  above,  the  importance  of  inter-impurity  interactions  became 
apparent  through  the  study  of  heavy  fermion  compounds,  where  the  impurity 
separation  R is  less  than  the  radius  of  the  cloud  of  localized  electrons  which  screen 
each  impurity  (screening  length  £),  R < £ ~ vp/ksTK-  these  compounds 
the  overlap  of  the  two  clouds  of  localized  electrons  formed  around  two  separate 
impurity  sites  causes  an  effective  coupling  between  the  impurities,  known  as  the 
Ruderman-Kittel-Kasuya-Yosida  (RKKY)  interaction  [54,  33,  74].  The  presence 
of  an  inter-impurity  interaction  can  interfere  with  the  ability  of  the  Kondo  effect 
to  screen  the  impurities.  The  two-impurity  Kondo  model  has  been  used  to  study 
the  interplay  of  these  two  processes  [31,  32], 

The  two-impurity  Kondo  model  describes  two  local  moments  coupled  to  a 
noninteracting  conduction  band, 

2 

H2imP  = Y1  eki ack,a  - J ^2  Sj  • s(r,).  (1-18) 

k<r  j= 1 

As  in  the  one-impurity  model  the  conduction  band  operators  c k<j  annihilate  Bloch 
states  of  energy  e k with  momentum  k and  spin  z-component  a. 

The  exchange  coupling  between  the  two  impurities  and  the  conduction  elec- 
trons produces  an  RKKY  interaction 

Hrkky  Rkky^l  ' S2,  (1.19) 

through  the  virtual  exchange  of  conduction  electrons  between  the  impurity  sites. 
As  with  the  exchange  coupling  parameter  a positive(negative)  RKKY  coupling 
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represents  a(n)  (anti)ferromagnetic  interaction.  The  RKKY  coupling  depends  both 
on  the  impurity  separation  R and  the  exchange  coupling  parameter  J: 

Rkky  J2  cos(kFR)/(kFR)3,  (1-20) 

where  kp  represents  the  Fermi  wavevector. 

Numerical  renormalization-group  calculations  find  the  strength  of  the  RKKY 
coupling  relative  to  the  Kondo  temperature  determines  the  ground-state  impurity- 
spin  configuration  at  low  temperatures  [31,  32].  For  ferromagnetic  or  weak  antifer- 
romagnetic RKKY  couplings,  the  Kondo  effect  screens  the  impurity  spins  before 
the  RKKY  interaction  can  correlate  them.  On  the  other  hand,  when  the  RKKY 
interaction  is  strongly  antiferromagnetic  the  impurity  spins  form  a spin-singlet, 
which  decouples  from  the  conduction  band  before  the  Kondo  effect  can  screen  the 
impurities.  In  both  of  these  regimes  the  low-energy  excitations  of  the  two-impurity 
model  are  described  by  a Fermi  liquid. 

Under  conditions  of  strict  particle-hole  symmetry,  there  is  an  abrupt  transi- 
tion between  the  two  regimes  described  above  at  a universal  ratio  of  the  RKKY 
coupling  to  the  Kondo  temperature  ( Irkky/Tp  % ! — 2.4  independent  of  the  value  of 
J).  At  this  transition  neither  the  Kondo  effect  nor  the  RKKY  interaction  is  able 
to  dominate  the  other  and  the  system  has  residual  interactions  down  to  T = 0. 
The  resulting  non- Fermi- liquid  properties  are  closely  analogous  to  those  of  the 
two-channel  Kondo  model  (see  below). 

Away  from  particle-hole  symmetry  there  is  no  abrupt  transition  and  the 
low-temperature  behavior  is  always  that  of  a Fermi  liquid. 
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Table  1.1:  Some  materials  which  display  non- Fermi  liquid  properties  at  low  tem- 
peratures and  the  mechanisms  proposed  to  explain  these  properties. 


Alloy 

Proposed  Mechanism 

Reference 

Thi_xUxRu2Si2 

Two- channel 

[5,  6] 

Thi_xUxPd2Si2 

Two-channel 

[7] 

Lcti —x  Cgx  Cu2Si2 

Two- channel 

[13] 

UCu5_xPdx 

Disorder 

[14,  20] 

CeCu6-xAux 

QCP 

[40] 

CePd2Si2 

QCP 

[25] 

1.4  Non-Fermi  Liquids 

Table  1.1  lists  a sample  of  /-moment  systems  which  generate  non- Fermi- 
liquid  behavior  at  low  temperatures.  Three  mechanisms  have  been  proposed  to 
describe  these  non-Fermi-liquid  systems:  multichannel  and  disordered  Kondo  ef- 
fects as  well  as  proximity  to  a quantum  critical  point.  Below  we  outline  each  of 
these  mechanisms.  We  also  propose  magnetic  frustration  as  a fourth  route  to  non- 
Fermi-liquid  behavior  and  introduce  a simple  model  which  incorporates  frustration 
to  test  the  feasibility  of  this  hypothesis. 

1.4.1  Multichannel  Kondo  Models 

Multichannel  Kondo  models  have  the  same  form  as  the  Kondo  models  pre- 
sented above,  except  that  the  local  moments  couple  to  more  than  one  conduction 
band.  The  Hamiltonian  which  describes  a one-impurity  M-channel  Kondo  model 


has  the  form 
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M 

Hmulti  ~ 52  (52  Ck<rmCk<rm  — J S • Sm).  (1-21) 

m—  1 kcr 

Here  sm  measures  the  spin  of  the  conduction  electrons  in  channel  m at  the  impurity 
site  (taken  to  be  the  origin),  i.e., 

Sm  = SE  clamaayCkc'm.  (1.22) 

kk'  at t' 

Multichannel  Kondo  models  were  introduced  by  Nozieres  and  Blandin  [49], 
who  showed  the  low-temperature  behavior  of  these  models  depends  primarily  on 
the  size  of  the  impurity  spin  relative  to  the  number  of  channels.  These  authors 
divided  multichannel  models  into  three  categories:  underscreened,  screened,  and 
overscreened.  In  an  underscreened  model  (M  < 2 5)  the  number  of  channels  is 
insufficient  to  screen  the  impurity  spin.  A residual  local  moment  with  spin  S—M/2 
remains  after  the  exchange  coupling  goes  to  strong  coupling.  In  a screened  model 
(M  = 25),  the  Kondo  effect  completely  quenches  the  impurity  spin.  The  original 
Kondo  model  (M  = 1,  5 = 1/2)  typifies  screened  models.  Both  underscreened  and 
screened  models  generate  Fermi-liquid  behavior  at  low  temperatures. 

By  contrast,  overscreened  models  (M  > 25)  can  exhibit  non-Fermi-liquid 
properties  at  low  temperatures.  At  strong  coupling,  a residual  local  moment  with 
spin  M/2-5  forms  from  the  overcompensated  impurity  spin  and  M electrons 
which  couple  to  the  impurity.  Nozieres  and  Blandin  showed  this  residual  local 
moment  has  a weak  antiferromagnetic  coupling  to  the  conduction  band.  Since  a 
weak  antiferromagnetic  coupling  is  unstable  and  scales  to  stronger  coupling  (just 
as  in  the  conventional  Kondo  model),  the  strong-coupling  limit  must  be  unstable. 
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By  the  same  token,  the  weak-coupling  configuration  of  an  unscreened  impurity 
must  also  be  unstable.  Therefore,  Nozieres  and  Blandin  concluded  the  exchange 
coupling  between  the  unscreened  impurity  and  the  conduction  band  must  take  an 
intermediate  value. 

The  two-channel  S = 1/2  Kondo  model  is  the  most  widely  studied  over- 
screened model.  Possible  applications  of  this  model  include  the  quadrupolar  Kondo 
effect  due  to  U4+  ions,  for  example  in  Uo.gTho.iBe^  [4],  and  magnetic  two-channel 
materials  formed  from  Ce3+  ions,  for  example  La1_ICexCu2Si2  [13,  12]. 

The  low-temperature  properties  of  overscreened  multichannel  models  have 
been  studied  using  a variety  of  techniques.  We  focus  on  the  results  derived  through 
Bethe  ansatz  [15,  55,  68]  and  conformal  field  theory  calculations  [41]  summarized 
by  Cox  and  Zawadowski  [19].  The  susceptibility,  linear  specific-heat  coefficient, 
and  Wilson  ratio  of  the  overscreened  S — \ M-channel  Kondo  models  are  given 
by 


Ximp  OC  {T/TkY  \ 

limp  oc  (T/TKy-\  (1.23) 

(2  + 2/M)(2  + M)2 


Here  r = 4/(2  + M). 

The  two-channel  model  (M  = 2)  represents  a special  case  with  a logarithmic 
low-temperature  susceptibility,  a logarithmic  linear  coefficient  of  the  specific  heat, 
and  a Wilson  ratio  of  8/3.  The  presence  of  an  unquenched  object  formed  by 
the  impurity  and  conduction  electrons  at  low  temperatures  leads  to  a very  large 
susceptibility. 
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It  should  be  noted  non-Fermi-liquid  behavior  only  occurs  in  multichannel 
Kondo  models  when  each  channel  experiences  the  same  exchange  coupling  with  the 
impurity,  Jm  = J.  Asymmetry  in  these  couplings  tips  the  competition  between  the 
channels  of  electrons  trying  to  screen  the  impurity  in  favor  of  the  channels  with  the 
largest  (in  magnitude)  exchange  and  generally  leads  to  conventional  Fermi-liquid 
physics. 

1.4.2  Kondo  Disorder 

Disorder  in  systems  containing  magnetic  impurities  leads  to  a distribution  of 
Kondo  temperatures.  If  the  distribution  of  TK' s is  sufficiently  broad,  unquenched 
local  moments  will  survive  to  arbitrarily  low  temperatures.  It  has  been  shown  [44], 
in  this  case  the  system  can  behave  as  a non-Fermi  liquid. 

One  can  introduce  disorder  into  a Kondo  model  by  varying  the  exchange  cou- 
pling randomly  from  site  to  site.  The  exact  form  of  the  distribution  of  exchange 
coupling  parameters  depends  on  the  low-temperature  properties  one  wishes  to 
reproduce.  The  disordered  Kondo  model  has  been  solved  using  dynamical  mean- 
field  theory  analysis  (which  neglects  any  effect  of  RKKY  interactions)  [44].  These 
calculations  reproduce  the  linear  dependence  on  temperature  of  the  resistivity  ob- 
served in  UCu5_xPdx  as  well  as  the  logarithmic  dependence  on  temperature  of  the 
susceptibility  and  the  linear  coefficient  of  the  specific  heat. 

1.4.3  Quantum  Critical  Points 

Proximity  to  a quantum  critical  point  provides  a third  mechanism  generating 
non-Fermi-liquid  behavior  in  heavy  fermion  compounds.  A quantum  critical  point 
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marks  a T = 0 phase  transition  produced  by  the  divergence  of  long-wavelength 
fluctuations  in  some  order  parameter.  In  /-moment  systems  such  transitions  may 
be  induced  by  the  application  of  pressure  or  changes  in  chemical  composition.  For 
example,  CeCue-xAu*  orders  antiferromagnetically  only  below  the  Neel  tempera- 
ture T/v-  Under  ambient  pressures  it  is  found  that  T/v  — ■>  0 as  x — > xc  ^ 0.1  [40]. 
For  x close  to  xc  the  Neel  temperature  can  also  be  driven  to  zero  through  the 
application  of  pressure.  The  correlations  induced  by  long-range  fluctuations  in  the 
magnetic  order  parameter  produce  non-Fermi-liquid  behavior  in  a finite  tempera- 
ture (T  < 2.5K)  phase  which  separates  the  antiferromagnetic  and  paramagnetic 
phases. 

1.4.4  Magnetic  Frustration 

Magnetic  frustration  is  a feature  of  many  heavy  fermion  systems  [38,  46]. 
It  forms  a fourth  possible  mechanism  for  non- Fermi-liquid  behavior  in  materials 
which  incorporate  interactions  among  the  local  moments.  Magnetic  frustration 
occurs  when  interactions  among  the  spin  degrees  of  freedom  of  a system  prevent 
the  individual  local  moments  from  each  accessing  their  ground  state. 

A toy  model  of  magnetic  frustration  is  provided  by  three  spin- 1 Heisen- 
berg spins  arranged  to  form  an  equilateral  triangle  (see  Fig.  1.1)  and  coupled 
antiferromagnetically.  One  can  describe  the  interactions  in  this  model  using  the 
Hamiltonian 


Hfms  = -IoY,Si ' SJ’ 

i<j 


(1.24) 
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Figure  1.1:  Spatial  configuration  of  three  spin-|  impurities,  each  constrained  to 
point  either  up  or  down.  Note  that  it  is  impossible  for  each  spin  to  point  in 
the  opposite  direction  to  both  of  its  neighbors.  This  is  an  example  of  “magnetic 
frustration.” 

where  i and  j label  the  impurity  sites  and  Iq  < 0.  In  this  example  it  is  favorable 
for  each  spin  to  line  up  antiparallel  to  its  two  neighbors.  As  it  is  impossible  for  all 
spins  to  achieve  this  state,  the  lowest-energy  configuration  accessible  to  the  model 
includes  one  frustrated  bond  where  one  pair  of  spins  remain  parallel.  The  model 
can  relieve  the  frustrated  bond  by  flipping  a spin  and  accessing  one  of  the  other 
spin- 1 configurations.  However  this  will  create  a new  pair  of  parallel  spins. 

Candidate  non-Fermi-liquid  materials  for  which  magnetic  frustration  might 
be  the  applicable  mechanism  incorporate  local  moments  which  form  a hexagonal 
or  Kagome  lattice  and  couple  antiferromagnetically.  Examples  of  this  type  of 
system  include  UNi4B  [48]  and  CePd2Al3  [59].  In  these  materials  the  geometry  of 
the  impurity  configuration  leads  to  frustrated  impurity  spins  at  low  temperatures. 
A simple  model  of  this  effect  is  a Kondo  Hamiltonian  for  three  spins  arranged 
as  in  Fig.  1.1.  When  the  RKKY  interaction  is  strongly  antiferromagnetic,  four 
degenerate  impurity  states  of  total  spin  one-half  (illustrated  in  the  last  four  lines 
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of  Table  1.2)  couple  to  three  degenerate  channels  of  conduction  electrons  (s  waves 
centered  on  each  impurity  site).  This  is  a higher  level  of  degeneracy  than  is  present 
in  the  two-channel  Kondo  model  and  is  highly  suggestive  of  the  existence  of  non- 
Fermi-liquid  behavior. 

1.5  A Three-Impurity  Kondo  Model 

We  now  formulate  a Kondo  model  motivated  by  the  three-impurity  Heisen- 
berg model  presented  in  the  previous  section.  Our  model  describes  a noninter- 
acting conduction  band  coupled  to  three  spin-|  impurities  through  exchange  and 
potential  scattering  interactions: 

Hump  = E ekCkaCk<T  + E^fo)  “ JS1  ’ S(ri)]‘  (L25) 

k,cr  j=l 

We  take  the  conduction  band  to  be  isotropic  in  momentum  space  and  to  extend 
over  a range  in  energy  ±D  centered  on  the  Fermi  energy  (ep  — 0).  The  con- 
duction electron  operators  c kCT  annihilate  Bloch  wave  states  of  momentum  k with 
spin  z-component  a.  These  operators  satisfy  the  fermionic  anticommutation  rela- 
tions defined  for  the  one-impurity  model,  Eq.  (1.5).  The  electronic  spin  operator 
s(r)  measures  the  electronic  spin  while  n( r)  measures  the  number  of  conduction 
electrons  at  r.  The  operators  Sj  act  on  the  spin-|  impurity  at  site  j. 

We  place  the  impurities  at  the  vertices  of  an  equilateral  triangle  in  the  x-y 
plane  and  centered  on  the  origin,  Vj  = (R/ \/3,  2nj /3,  7t/2)  in  spherical  coordinates. 
In  this  case,  the  Hamiltonian  given  in  Eq.  (1.25)  generates  an  effective  RKKY  inter- 
action between  the  impurities  having  exactly  the  same  form  as  given  by  Eq.  (1.24), 
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Table  1.2:  A basis  for  the  spin  state  of  three  spin-^  impurities.  Notations  such 
as  | TIT)  show  the  orientation  of  each  individual  impurity  spin,  while  the  entire 
configuration  is  denoted  \hS^ot).  Here  h is  the  helicity  of  the  state  as  defined  in 
Section  2.1  and  St2ot  represents  the  total  2-component  of  the  spin  of  the  state,  with 
Tf  representing  S(2ot  = §,  T representing  Sztot  = j representing  Szot  = —5,  and  Jj. 
representing  S20(  = — §. 


state  h Stot  Slot 

|o  fi)  = ITTT)  0 I 1 

1°  T)  = 73  ( IITT)  + I TIT)  + ITU))  5 

|0|)  = 73  (ITU)  + UU)  + lilt))  -5 

|01T)  = | ill)  -1 

|1T>  = 7s(  IITT)  + ei27r/3l  TIT)  + e-i27r/3|TTD)  1 \ \ 

IU)  = 75(ITU)  + ei2ir/3|  TTI)  + e-^lUT))  -5 

l2T)  = ^(UTT)  + e_i27r,/3|  TIT)  + ei27r/3|  TTi) ) 2 5 2 

|2i)  = 7s(ITU)  + e'i27r/3UTl)  + ei27r/3|UT)) 
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where  the  sign  of  /o(=  Irkky)  oscillates  with  increasing  impurity  separation  R.  For 
antiferromagnetic  70,  magnetic  frustration  is  built  into  this  model. 

It  is  useful  to  note  the  similarities  and  differences  we  expect  to  find  between 
our  three-impurity  model  and  the  two-impurity  model  studied  previously  [31,  32]. 

The  Kondo  temperature  and  the  RKKY  coupling  represent  the  relevant  en- 
ergy scales  in  both  models.  As  in  the  two-impurity  case,  the  RKKY  interaction 
and  the  Kondo  effect  should  not  compete  in  the  ferromagnetic  RKKY  regime  of 
the  three-impurity  model.  When  the  RKKY  interaction  dominates,  the  impurities 
should  form  a spin-|  configuration  which  the  conduction  band  screens.  When  the 
Kondo  effect  dominates,  the  conduction  band  should  screen  the  individual  impurity 
spins.  In  either  case  the  low-temperature  properties  should  be  the  same. 

In  the  regime  of  strong  antiferromagnetic  RKKY  coupling,  we  do  not  expect 
the  Kondo  effect  and  RKKY  interaction  to  commute.  When  the  Kondo  effect 
dominates  the  conduction  electrons  should  screen  the  impurity  spins  separately. 
When  the  RKKY  interaction  dominates  the  impurities  should  form  frustrated  spin- 
1 configurations,  which  then  couple  to  the  conduction  electrons.  As  discussed 
above,  this  seems  likely  to  lead  to  non- Fermi-liquid  behavior. 

1.6  Outline 

The  remainder  of  this  dissertation  describes  the  results  of  our  investigation 
of  the  three-impurity  Kondo  model  introduced  in  Section  1.5. 

In  Chapter  2 we  derive  the  numerical  representation  of  the  model  which  we 
use  to  analyze  its  low-temperature  properties.  We  identify  the  symmetries  of  the 
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model  and  discuss  how  they  may  be  exploited  when  calculating  the  low-energy 
spectra  within  the  numerical  representation. 

In  Chapter  3 we  identify  the  fixed  points  which  determine  the  low-temperature 
physics  of  the  three-impurity  model  and  describe  the  phase  diagram  of  the  model. 

In  Chapter  4 we  describe  the  results  of  more  detailed  analysis  of  the  regimes 
identified  in  Chapter  3.  When  possible  we  calculate  the  low-temperature  properties 
of  these  regimes. 

We  conclude  with  a summary  of  our  results  and  a discussion  of  possible 
avenues  for  future  research. 


CHAPTER  2 

NUMERICAL  FORMULATION  OF  THE  PROBLEM 


In  this  chapter  we  describe  how  we  transform  the  three-impurity  Kondo 
model  introduced  in  Section  1.5  into  a form  suitable  for  numerical  analysis.  We 
describe  the  symmetries  of  this  representation.  We  also  outline  important  aspects 
of  calculating  the  low-energy  spectra  within  the  numerical  representation. 

2.1  The  Discrete  Representation 

In  this  section  we  explain  how  we  formulate  a discrete  representation  of  the 
three-impurity  model.  This  involves  three  stages:  (1)  The  model  is  mapped  from 
three  spatial  dimensions  to  an  effective  one-dimensional  representation.  (2)  The 
continuum  of  conduction-band  states  is  discretized.  (3)  The  discretized  model  is 
cast  into  tridiagonal  form  using  the  Lanczos  method. 

2.1.1  One-Dimensional  Hamiltonian 

One  can  map  the  three-impurity  Kondo  model  onto  an  effective  model,  hav- 
ing one  spatial  coordinate:  the  radial  distance  from  the  center  of  the  triangle 
formed  by  the  impurities.  This  mapping  converts  the  original  conduction  band 
into  three  degenerate  conduction  bands,  each  formed  from  conduction  electrons 
with  a distinct  rotational  symmetry  or  “helicity.”  A helicity-/*  state  acquires  a 
multiplicative  phase-factor  of  ei2nh /3  (h  — 0, 1,2)  when  the  coordinate  system  is 
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rotated  by  27t/3  radians  about  the  origin.  Details  of  the  transformation  to  one 
dimension  are  given  in  Appendix  A.  The  resulting  Hamiltonian  takes  the  form 

H/D  = I deJ2  ^C\haCeha  + 2 Qq  FhFh>  53  fohafoh'a' 

•'  h,<x  hh'  aa' 

x[V5h,h'5*,<r'  - J\dc T,<7'  • s (h'-h)]-  (2-1) 

The  electronic  operators  c£hc  annihilate  states  of  reduced  energy  £,  helicity  h,  and 
spin  z-component  a.  The  reduced  energy  e represents  a dimensionless  parameter- 
ization of  the  energy  of  a conduction  electron  with  momentum  k, 


£ = tk/D.  (2.2) 

The  parameters  V and  J parameterize  potential  scattering  and  exchange 
interactions  respectively,  while  go  is  the  density  of  states  at  the  Fermi  energy.  The 
impurity  spin  operators 


Sh  = Y,ei2nhj/%  (2.3) 

j 

act  on  the  three  spin-|  impurities  to  increase  the  helicity  of  their  collective  con- 
figuration by  h.  The  operators  foha  couple  the  conduction  band  and  impurity 
configurations  and  are  defined  by 


Fhfoha  — yj^lg  j ^ ^ 9 h(keR)cehai 


(2.4) 


where  the  gF s satisfy 
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9h(x)  = - 

|[1  - sin(x)/x]  to  =1,2. 


|[1  + 2sin(x)/x]  to  = 0, 


(2.5) 


By  construction  the  fohaS  satisfy  fermionic  anticommutation  relations, 


(2.6) 


These  relations  and  the  shape  of  the  conduction  band  define  the  form  factors 


is  fully  equivalent  to  the  model  defined  in  Section  1.5. 

2.1.2  Logarithmic  Discretization 

In  order  to  solve  the  one-impurity  Kondo  model  Wilson  divided  the  conduc- 
tion band  into  a sequence  of  energy  intervals  whose  range  decreases  exponentially 
as  a function  of  an  interval  index  m [72], 


(2.7) 


With  this  definition  F%  — » 1/3  as  R — > oo. 

We  have  made  no  approximations  in  deriving  this  Hamiltonian,  so  Eq.  (2.1) 


A-(m+i)  < |£|  < A-m_ 


(2.8) 


Here  the  interval  index  takes  all  nonnegative  integral  values  (m  > 0)  and  one 
defines  the  discretization  factor  A to  be  greater  than  one.  Wilson  showed  that 
by  replacing  the  continuum  of  conduction  band  states  by  a discrete  subset,  one 


28 


can  introduce  an  artificial  separation  of  energy  scales  into  the  “hopping”  coeffi- 
cients of  the  tridiagonal  form  of  the  one-impurity  Hamiltonian  (see  below).  This 
allows  a convergent  approximation  to  the  infinite-chain  problem  described  by  the 
tridiagonal  Hamiltonian  using  finite-length  chains  which  correctly  reproduces  the 
impurity  contribution  to  the  properties  of  tlje  model.  We  use  the  same  process  of 
logarithmic  discretization  in  the  three-impurity  model  and  study  the  properties  of 
the  finite-length  chains  which  result. 

We  define  a complete  set  of  destruction  operators  ampha  ( bmphc ) and  an  as- 
sociated set  of  orthonormal  functions  ipmPh(£)  which  vanish  for  any  e outside 
the  mth  positive(negative)  bin.  Here  functions  with  a +(— ) superscript  are  de- 
fined for  positive(negative)  e and  the  harmonic  index  p takes  all  integral  values 
(-00  < p < oo ).  Given  such  a basis,  one  can  write 


OO  OO 

ceha  — JZ  rnph(£)amphc  T t^mph  (e )^mpha]  ■ (2-9) 

m=0  p=— oo 

The  particle  ( ampha ) and  hole  (bmpha)  operators  with  helicity  h and  spin  z-component 
a are  superpositions  of  the  c£hc  operators  within  each  interval, 

amph<T  = A~\  f de{^  {£))*ceha  (2.10) 

J A-'m+n 

and 

= B-‘h  y_A_m  <fe«w(£))W  (2.11) 

These  particle-hole  operators  satisfy  fermionic  anticommutation  relations, 
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WmphtT^ln'p'h'cT'}  ~ {fynp/KTi  ^L'p'/iv)  — <5m,m'<5p,p'<5/i,/i'<Ws  (2-12) 

{®rapki  ^n'p'/i'ff'}  (2-13) 


The  key  step  in  solving  the  Kondo  model  is  to  choose  the  p — 0 function  so 


that  it  has  the  same  energy  dependence  as  wh{e)  = \Jg(e)/(2Qo)gh(keR): 


and 


where 


TpmOh(e)  = i 


TpnO  h(£)  = 


A~lwh{£)  for  A~(m+1)  < £ < A_m, 
0 otherwise; 


B-hV(e)  for  - A“m  < e < -A-(m+1), 
0 otherwise; 


(2.14) 


(2.15) 


A2 


mh 


r>2 
^ mh 


rA-m 

/ i ^dswh(£), 

I — A~(m+0 

/ d£w2h(e). 


(2.16) 

(2.17) 


With  these  definitions  the  on-site  electronic  operators  take  the  form 


foha  — ^ J^Omh^mOha  T ^Om/i^mO/w]' 


(2.18) 


To  ensure  orthonormality, 
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-A-mh 

V'Omh  nr  5 ^Om/i 

V 


Bmh 

V2Fh 


(2.19) 


where  Fh  is  defined  in  Eq.  (2.7).  This  leads  to  a discretized  representation  of  the 
conduction  band  Hamiltonian, 


Hc/D  — |(1  + A *)  \Qmhamhaarnhcr  + Pmh^rnha^rnhcpi  (2.20) 

nha 

with  the  coefficients 

94-2  r\-m 

amh  = 1 ,T- 1 / , +n  (2-21) 

1 + A 1 7A-(m+l) 

od- 2 

Pmh  = TV TT  / ew2h{e)d£.  (2.22) 

1 + A-1  7-A-m 

2.1.3  Tridiagonal  Hamiltonian 

We  now  use  the  Lanczos  method  [39]  to  transform  Hc  as  defined  in  Eq.  (2.20) 
into  a tridiagonal  or  “hopping”  form, 

Hc/D  — ^2  [£nhfnhafnh<r  + tnh(flh<rfn+lh<r  + fl+lhafnho)\,  (2.23) 

n,h,a 

in  which  electrons  hop  along  a linear  chain  of  sites  terminated  at  one  end  by  the 
impurity  site.  This  semi-infinite  chain  of  sites  forms  a tight-binding  representation 
of  the  logarithmically  discretized  conduction  band. 

The  / operator  which  represents  an  electron  at  site  n has  the  form 

OO 

fnha  — ^ “1“  ^nmhf^mha)  j (2.24) 

m= 0 
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where  the  coefficients  unmh  and  vnmh  define  an  orthogonal  transformation  from  the 
basis  of  the  discretized  conduction  band  to  the  tight-binding  basis. 

The  kinetic  energy  coefficients  enho,  tnho  and  the  coefficients  unmh,  vnmh  are 
determined  from  Hc  — 2//c/(l  + A-1)  and  the  specification  of  foha  in  Eq.  (2.18) 
via  a set  of  recursion  relations: 


£nh  = (fnha\Hc/D\fnha),  (2.25) 

4 = \[(HC/D  - enh)\ fnha)  - tn-lh\fn-ih*)}\2,  (2-26) 

with 

\fn+lhc)  = 4[(-^c/ D — £nh)\fnhc)  — tn-\h\ f n-\hcr)\-  (2.27) 

In  this  basis,  the  operators  fohc  combine  electrons  which  span  the  entire 
range  of  reduced  energy  |e|  < 1.  The  operators  on  site  n,  on  the  other  hand, 
combine  electrons  confined  to  a range  of  reduced  energy  |e|  < A-n//2.  As  a result 
the  states  annihilated  by  the  /o/ia’s  are  the  most  localized  about  the  impurity  site, 
with  the  states  annihilated  by  the  fnha' s becoming  less  localized  with  increasing  n. 

One  can  formulate  the  discretized  Hamiltonian  in  terms  of  a sequence  of 
Hamiltonians  H ^ acting  on  finite-length  chains  of  sites, 


H/D  = (Hc  + Hint)/D  = lim  |(1  + A-^A"^-1^2^.  (2.28) 

N — KX) 


A recursion  relation  defines  i7/v+i> 
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Hn+i  — A1(/“//iv  + y~l[^jv+i/i/jv+ifcg/jv+ifc<7 


h,cr 


+^Nh{fNhafN+lha  + flr+lkaf^ha)]- 


(2.29) 


Here 


and 


TT  Hint 

(1  + A-1)  D 
tNh  = A ^N~l^2£Nh, 

t,Nh  = hN^2tNh- 


(2.30) 


(2.31) 

(2.32) 


Thus,  the  lowest  energy  scale  of  interactions  represented  by  the  Hn+i  [given  by 

CN+lhfN+ihtrfN+lha  + £ Nh(flrhafN+lh<r  + fN+lhafNha)]  IS  of  Order  1. 

This  sequence  of  Hamiltonians  reproduces  the  Hamiltonian  acting  on  the 
continuous  one-dimensional  conduction  band  in  the  continuum  limit  (i.e.,  A — * 1), 
as  the  iteration  number  N goes  to  infinity. 


2.1.4  Independent  Couplings  Approximation 

We  use  an  independent  couplings  approximation  [31]  to  reduce  the  dimension 
of  the  parameter  space  which  we  explore  numerically.  The  discretized  Hamiltonian 
defined  in  Eq.  (2.28)  contains  an  infinite  number  of  parameters:  qqJ,  QqV , Fh, 
and  enh . The  overwhelming  majority  of  these  parameters  are  associated  with  terms 
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in  the  Hamiltonian  which  are  irrelevant  in  the  renormalization-group  sense.  That 
is,  the  strength  of  these  couplings  decreases  as  a power  of  the  temperature  under 
renormalization  and  vanishes  as  T — > 0.  The  only  exceptions  are  the  potential 
scattering  and  exchange  terms  entering  Ho,  which  are  marginal  and  marginally 
relevant,  respectively.  Marginal  couplings  remain  constant  under  renormalization, 
while  marginally  relevant  couplings  grow  logarithmically  with  temperature.  These 
couplings  determine  the  low-energy  properties  of  the  three-impurity  Kondo  model. 
The  sole  effect  of  the  irrelevant  interactions  is  to  shift  the  effective  value  of  marginal 
couplings  and  possibly  to  break  certain  symmetries  exhibited  among  the  marginal 
couplings. 

The  full  discretized  Hamiltonian,  Eq.  (2.28),  contains  just  five  independent 
marginal  or  marginally  relevant  couplings:  qqJFq,  QqJ F± , Q0JF1F2,  QqVFq,  and 
QoVF?.  However,  one  can  show  using  an  extension  of  poor  man’s  scaling  [9]  for  the 
three-impurity  model  that  under  renormalization  irrelevant  interactions  can  split 
the  terms  multiplied  by  eoJF?  into  two  independent  groups,  each  of  which  has  its 
own  effective  coupling.  We  are  therefore  left  with  a minimal  set  of  six  marginal 
parameters  describing  the  model.  The  preceding  arguments  indicate  that  one  can 
explore  the  entire  range  of  behaviors  of  the  three-impurity  model,  including  the 
effect  of  all  possible  band  structures  g(e)  and  impurity  separations  R,  using  an 
independent  couplings  approximation  which  fixes  the  kinetic  energy  coefficients 
but  allows  the  six  marginal  couplings  to  take  all  possible  values. 

We  construct  Hc  in  the  independent  couplings  approximation  from  kinetic 
energy  coefficients  describing  a flat,  symmetric  conduction  band  and  infinite  im- 


34 


purity  separation  [p(f)  = Qo  and  R = oo],  in  which  limit  enh  and  £nh  reduce  to  the 
coefficients  for  the  one  impurity  model  [72], 


^nh 


6n  0} 


1 - 


(2.33) 


inh  —>  fn  = ^ _ A_(2„+l)^1  _ (2n+3))]l/2' 

The  interactions  between  the  impurities  and  conduction  electrons  have  es- 
sentially the  same  form  as  those  described  by  the  full  discretized  Hamiltonian: 


Here 


Hint/D  = [JooSoo  + JnSn  + J11822]  • So 

+ [<7oiSoi  + </ois20  + <7l2s12]  ■ Si 
+ ['7ois02  + ^OlSlO  + ^12s2l]  ' S2 
+ J2V^h-  (2-34) 

h,<r 


and 


Shh'  — foha  2 foh'i 


a ,cr' 


(2.35) 


nh  = Y,  fohJohc-  (2-36) 

O 

Here  the  form  factors,  density  of  states,  and  exchange  coupling  are  absorbed  into 
one  of  four  couplings: 
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—QoJFhFh' 


Joo  for  h = h'  = 0 
Jn  for  h = ti  £ { 1,2} 

J12  for  h ^ h1  e {1,  2} 

Joi  for  h = 0,  h!  = 1, 2 or  h'  = 0,  h = 1, 2. 


(2.37) 


Note  that  the  couplings  representing  QqJF?  and  goJF^  are  identical  due  to  a 
symmetry  of  the  three-impurity  model,  whereas  the  coupling  representing  Q0JF1F2, 
although  identical  to  these  couplings  initially,  can  assume  different  values  under 
renormalization  and  must  therefore  be  treated  separately. 

We  normalize  the  values  of  the  Ji  s to  ensure  Eq.  (2.34)  correctly  reproduces 
Eq.  (2.1)  in  the  limit  R — ► 00: 


We  construct  the  coefficients  of  the  allowed  potential  scattering  processes  from  the 
form  factors  and  the  coupling  parameter  V, 


These  coefficients  also  determine  the  higher  order  interactions  represented 


Joo  + 2 Jn  + 2J^2  + 4Jqi  — (2qqJ)  . 


(2.38) 


QoVFl  — 


Vo  for  h = 0 


(2.39) 


Vi  for  h = 1,2. 


by  the  independent  couplings  Hamiltonian.  For  example,  the  effective  RKKY 
coupling  is 


(2.40) 
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As  we  discuss  in  Section  1.3,  the  strength  of  the  RKKY  parameter  depends  on 
impurity  separation.  Therefore,  Eq.  (2.40)  demonstrates  that  varying  the  exchange 
couplings  varies  the  effective  impurity  separation  represented  by  Hint. 


In  this  section  we  present  the  issues  relevant  to  the  process  of  diagonalizing 
our  numerical  representation  of  the  three-impurity  model.  We  first  introduce  the 
symmetries  and  conserved  quantities  of  H jv-  We  then  describe  how  we  expedite 
the  process  of  evaluating  the  matrix  elements  of  H w by  using  the  conserved  quan- 
tities to  decrease  the  number  of  matrix  elements  of  Hn  which  need  be  evaluated. 
We  discuss  how  the  conserved  quantities  may  also  be  exploited  to  expedite  the 
diagonalization  process.  We  conclude  with  an  overview  of  the  errors  associated 
with  determining  the  eigenvalues  of  H pj. 

2.2.1  Symmetries  and  Conserved  Quantities 

Four  quantities  are  conserved  due  to  the  symmetries  of  the  numerical  repre- 
sentation of  the  three-impurity  model: 


2.2  Diagonalization  of  H ^ 


total  spin 


S = Sc  + So 


total  z-component  of  spin  Sz  = Sz  + Sq 


charge 

helicity 


<3 

h = [hc  + himp)  mod  3. 
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Table  2.1:  Impurity  spin  operators  for  the  independent  couplings  Hamiltonian.  See 
Table  1.2  for  an  explanation  of  the  notation  for  the  impurity  states. 

Sq  = | ( |o  t)(o  fr|  - |o  4)<o  ) + 5 ( |o  t)(o  1 1 — |o  i)(o  i 1 
+|1T)(1T|-|U)(U|  + |2T>(2T|-|2  1X2  11), 

5+  = v/3(|0lf)(0T|  + |0|)(0^|)  + 2|0T)(0I|  + |lT)(Ul  + |2T)(21|, 

Sq  = Vz  ( |0  T)(0  -ffl  + |0  41) (0  1 1 ) -h  2|0  D(0  r I + |1  0(1  1 1 + |2  1)(2  T I, 
s?  = |2 1)<1  II  - 12  T)(1  Tl  - 10 1X2 II  - 10  T)(2  T|  - II 1X01I-I1TX0TI, 

St  = v/3(  |0  1T)(2  Tl  - |1  i)<0  411 ) - 2|2  T)(l  il  + |0  r>(2  il  - |1  T)(0  il, 

Sf  = >/3(|l  T)(01H  - |0  41)<2  1 1 ) — 2|2  D(1  Tl  + |1  1)(0  Tl  - |0  1>(2  T|> 

= II 1X2  11- 11  TX2TI- 12 1X0 II  - 12  T)(0  Tl  - 10 1)(1 11-10  T)(1  Tl, 
St  = >/3  ( |0  Tr)(l  T I — 12  1) (o  4M ) — 2|1  T>(2  11  + |0  T)(l  II  — |2  T > <0  11, 

5;  = V3(  |2  T)(0  4H  - |0  41><1  Tl ) - 2|1 1X2  Tl-b  |2  D(0  Tl  - 10  l)(i  T|. 


The  operators  which  measure  these  quantities  for  an  arbitrary  eigenstate  generate 
the  symmetries  of  H^. 

We  define  the  electronic  spin  operators  Sc  in  terms  of  the  Pauli  matrices  used 
to  describe  the  exchange  interaction, 

OO 

Sc=  £ Y.fnhc\°°,°'fnho'-  (2-41) 

n=0,h  a, a' 

We  present  the  definitions  of  the  impurity  spin  operators  S in  Table  2.1. 

The  operator  which  measures  the  charge  of  an  eigenstate  of  takes  the 


form 
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Q=Y.  ~ 1/2).  (2-42) 

n= 0 h,ar 


In  the  absence  of  potential  scattering,  HN  is  also  symmetric  under  the  particle-hole 
transformations 


/„»„  — * (-1  r/L,,  (2-43) 

S/i  ‘ — > — (2-44) 

or  under  the  transformation 

u,  — (-ir(-i)1/J-'/V*>»-  (2-45) 

Here  a represents  the  spin  antiparallel  to  a.  We  use  this  particle- hole  symmetry  to 
formulate  charge  conservation  in  terms  of  “axial  charge”  or  “isospin”  [32].  Axial 
charge  I is  a vector  operator  whose  ^-component  maps  onto  the  charge  Q of  the 
model  [given  in  Eq.  (2.42)], 


Iz  = \Q.  (2.46) 

We  represent  the  transverse  components  of  axial  charge  by  raising  and  lowering 
operators  ^ — Ix  ± ily , where 


/+  = E E(-i)"At Aj.  r = (/+)'-  <2-47> 

n= 0 h 

Angular  momentum  commutation  relations  hold  for  the  components  of  axial  charge. 
The  operator  which  measures  the  electronic  helicity  takes  the  form 
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he  = Hr  Y,  hfhjnia  mod  3, 


(2.48) 


\n= 0 h,a 


while  the  helicity  of  the  impurity  configuration  is  given  by 


himp  = mod  3.  (2.49) 

Hn  is  also  symmetric  under  the  exchange  of  the  helicity- 1 and  -2  labels  of  all 
electronic  and  impurity  operators.  This  symmetry  results  in  a degeneracy  between 
pairs  of  helicity- 1 and  -2  eigenstates  whose  quantum  numbers  are  otherwise  the 
same. 

In  diagonalizing  H # we  label  all  the  eigenstates  of  by  their  total  spin, 
total  ^-component  of  spin,  2-component  of  axial  charge,  and  helicity.  In  the  absence 
of  potential  scattering  we  also  use  total  axial  charge  to  label  the  eigenstates. 


2.2.2  Numerical  Treatment 

We  base  our  solution  of  H n on  work  done  on  the  two-impurity  model  [31,  32]. 
Once  one  takes  into  account  the  different  symmetries  of  the  two-  and  three-impurity 
models,  one  uses  essentially  the  same  algorithm  to  analyze  these  models.  In  fact, 
we  use  the  same  computer  code  to  study  one-,  two-,  and  three-impurity  models. 

We  outline  in  Appendix  B how  we  exploit  spin  and  overall  charge  conserva- 
tion to  reduce  the  size  of  the  basis  of  H ^ and  how  we  simplify  the  evaluation  of 
matrix  elements  of  H^. 

Using  the  conserved  quantities  defined  in  the  previous  section,  we  can  break 
Hn  down  into  block-diagonal  form,  where  each  block  or  submatrix  of  Wv  is  formed 
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from  the  matrix  elements  of  states  with  the  same  total  spin,  helicity,  and  axial 
charge  when  particle-hole  symmetry  is  obeyed,  and  otherwise  from  the  same  total 
spin,  helicity,  and  2-component  of  axial  charge.  The  time  required  to  diagonalize 
Hn  as  a whole  depends  linearly  on  the  number  of  submatrices  of  HN,  while  the  time 
required  to  diagonalize  a submatrix  goes  as  the  cube  of  its  dimension.  Thus  the 
smaller  one  can  make  each  of  the  submatrices,  by  working  with  a greater  number 
of  conserved  quantities,  the  faster  the  process  of  diagonalization. 

Errors  in  calculating  the  eigenvalues  of  Hn  have  two  sources:  truncation 
of  the  basis  and  round-off  errors  introduced  during  numerical  diagonalization  of 
the  Hamiltonian.  Truncation  is  unavoidable,  because  the  full  basis  of  the  discrete 
Hamiltonian  grows  by  a factor  of  26  at  each  iteration.  After  only  a few  iterations  it 
becomes  impossible  to  carry  forward  every  eigenstate  of  iteration  N to  form  basis 
states  at  iteration  N + 1.  Instead  one  retains  the  M states  of  lowest  energy,  where 
M is  typically  in  between  500  and  2000. 

The  error  introduced  by  truncating  the  basis  of  H n tends  to  be  of  greater 
impact  than  the  round-off  error.  Once  the  maximum  number  of  retained  states  has 
been  reached  (around  iteration  5 typically),  one  keeps  only  a fraction  1/26  ~ 1.5% 
of  the  eigenstates  calculated  at  each  iteration.  This  is  a much  smaller  retained 
fraction  of  the  eigenstates  than  is  possible  for  the  single-impurity  model,  where 
the  basis  grows  only  by  a factor  of  22  at  each  iteration. 

Another  measure  of  the  effect  of  truncation  is  the  range  of  energies  of  the 
eigenstates  of  Hn  (measured  relative  to  the  ground-state  energy).  Recall  that 
due  to  the  rescaling  by  A1/2  at  each  iteration  of  Eq.  (2.29),  the  smallest  energy 
scale  in  Hn  is  of  order  unity.  In  the  single-impurity  Kondo  model  keeping  1000 
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states  ensures  an  energy  range  of  at  least  20.  By  contrast,  keeping  4000  states  in 
the  three-impurity  model  (a  number  much  higher  than  is  generally  practical  due 
to  the  enormous  computer  resources  required)  generates  an  energy  range  of  only 
10  or  so.  As  we  will  see  in  Chapter  4,  this  decrease  in  width  of  the  spectrum 
is  a major  obstacle  to  the  calculation  of  thermodynamic  properties  of  the  three- 
impurity  model. 

Round-off  error  tends  to  be  a smaller  source  of  error  than  truncation  error. 
However,  round-off  error  can  break  the  symmetry  between  states  which  should  be 
degenerate,  in  particular  ones  which  differ  only  in  the  interchange  of  the  helicity-1 
and  -2  label.  This  error  compounds,  so  after  40  to  50  iterations  these  splittings 
can  make  identification  of  fixed-point  symmetries  impossible. 

The  greatest  source  of  error,  however,  comes  from  logarithmic  discretization. 
Discretization  introduces  errors  in  the  results  calculated  for  the  continuum  limit 
from  the  energy  spectrum  generated  by  Hy.  We  compensate  for  this  error  by 
varying  the  value  of  A.  This  allows  us  to  estimate  the  errors  introduced  by  dis- 
cretization. In  general,  we  set  A equal  to  3.0.  This  seems  to  provide  the  best 
compromise  between  computational  efficiency  and  faithful  approximation  of  the 
continuum  limit. 


CHAPTER  3 

PRELIMINARY  RESULTS 

In  this  chapter  we  describe  the  different  low-temperature  regimes  of  our  three- 
impurity  Kondo  model  which  we  identify  from  the  energy  spectra  generated  by  the 
Hamiltonian  HN  defined  in  Eqs.  (2.28)-(2.32)  and  (2.34).  We  first  review  some 
pertinent  concepts  from  renormalization-group  theory,  in  particular  the  role  of  fixed 
points.  We  then  describe  the  spectra  and  stability  with  respect  to  perturbations  of 
each  fixed  point  of  the  three-impurity  model,  as  well  as  our  reasons  for  classifying 
each  fixed  point  as  Fermi-liquid  or  non-Fermi-liquid.  We  conclude  by  discussing 
the  full  phase  diagram  of  the  particle-hole-symmetric  model. 

3.1  Fixed  Points  and  Flows 

In  this  section  we  describe  how  we  discern  different  low-temperature  phases 
in  our  model  from  the  spectra  generated  by  the  discretized  Hamiltonian  derived 
in  the  previous  chapter.  We  start  with  a brief  discussion  of  renormalization-group 
transformations.  We  define  the  concepts  of  fixed  points  and  flows  generated  by  a 
renormalization-group  transformation.  We  also  establish  a relation  between  tem- 
perature and  the  iteration  number  N.  We  conclude  with  an  overview  of  the  phase 
diagram  of  the  three-impurity  model. 

A renormalization-group  (RG)  transformation  describes  a process  by  which 
one  calculates  the  properties  of  models  in  which  infrared  or  ultraviolet  divergences 
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occur.  The  RG  transformation  developed  by  K.  G.  Wilson  for  the  Kondo  prob- 
lem [73,  72]  provides  a nonperturbative  method  for  calculating  the  low-temperature 
spectra  of  the  model  by  incorporating  the  effects  of  interactions  on  different  energy 
scales. 

Following  Wilson,  we  base  the  RG  transformation  for  our  model  on  the  re- 
cursion relation  in  Eq.  (2.29),  with  the  ground-state  energy  subtracted  off: 


Here  E^+1  represents  the  lowest  eigenvalue  of  i/jv+i-  Formally,  we  can  define  a 
numerical  RG  transformation  RA  which  maps  the  eigenvalues  and  eigenstates  of 
Hn  onto  those  of  Hn+ i,  or  symbolically 


At  each  iteration  of  Eq.  (3.2),  the  eigenspectrum  of  iteration  N is  expanded 
by  a factor  A 1/2  (A  > 1),  then  new  terms  of  order  unity  are  added  to  form  i/jv+i- 
As  a result,  diagonalization  of  H^+i  introduces  splittings  into  the  spectrum  gen- 
erated at  iteration  N and  the  energy  levels  of  successive  spectra  bifurcate  as  one 
incorporates  more  degrees  of  freedom. 

3.1.1  Renormalization-Group  Concepts 


'0 

'N+ 1- 


(3.1) 


Hn+ i = RaHm. 


(3-2) 


One  of  the  most  important  aspects  of  the  numerical  solutions  is  the  existence 
of  fixed  points,  i.e.,  of  spectra  which  are  invariant  under  an  RG  transformation. 


44 


For  technical  reasons  discussed  below,  the  fixed  points  we  study  are  those  of  R A 
rather  than  R\,  i.e.,  we  study  fixed-point  spectra  of  Hamiltonians  H where 

R\H*  = H\  (3-3) 

In  practice,  identification  of  fixed  points  involves  iterating  the  RG  transforma- 
tion many  times  (typically  40-60)  until  the  eigenvalues  of  HN  agree  to  within  a 
prescribed  numerical  precision  with  those  of  Hyv-2- 

As  we  discuss  in  the  derivation  of  Hn,  the  operators  f^ha  annihilate  states 
which  span  a conduction  band  of  width  2A  D , where  the  on-site  electrons  fohcr 
span  a conduction  band  of  width  2D.  One  interprets  this  decrease  in  the  effective 
width  of  the  conduction  band  as  reflecting  a decrease  in  the  effective  temperature, 
where  the  thermal  energy  kBT  defines  the  upper  bound  on  excitations  energies. 
On  this  basis,  one  associates  each  iteration  with  a characteristic  temperature 


kBT  ~ A ~N/2D.  (3.4) 

Therefore,  we  see  that  fixed-point  spectra,  corresponding  to  the  limw—oo  Hn,  de- 
scribe the  zero-temperature  behavior  of  the  model  under  study. 

The  evolution  of  the  spectrum  of  HN  with  increasing  N can  be  visualized 
by  plotting  the  energy  of  low-lying  states  versus  iteration  number.  Since  odd  and 
even  iterations  tend  to  yield  very  different  spectra,  one  generally  only  plots  either 
odd- A or  even -N  spectra  on  the  same  “flow  diagram.”  Figure  3.1  shows  two 
regions  where  the  energy  levels  do  not  vary  appreciably  for  several  iterations.  The 
set  of  scale-invariant  levels  starting  around  iteration  46  represent  the  stable  fixed 
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Energy  Flows 


Figure  3.1:  Sample  flow  diagram  representing  the  even -N  energy  spectra  generated 
by  the  numerical  representation  for  A = 3.0,  Joo  = J\\  = 0.238943,  J n = 0.2,  and 
Joi  = 1-0.  This  diagram  presents  two  fixed  points:  an  unstable  fixed  point  from 
iteration  16  to  30  and  a stable  one  which  starts  at  iteration  46. 

point  describing  the  T = 0 state  of  the  system  (N  ->  oo).  The  intermediate  set 
of  scale-invariant  energy  levels,  generated  between  iteration  16  and  30,  identify 
an  unstable  fixed  point  near  which  the  initial  Hamiltonian  Ho  starts.  One  says 
the  initial  Hamiltonian  flows  under  renormalization  from  the  unstable  fixed  point 
toward  the  stable  fixed  point. 

The  designation  of  a fixed  point  as  stable  or  unstable  identifies  the  effect 
of  renormalization  on  perturbations  about  the  fixed-point  Hamiltonian  H* . Per- 
turbations about  stable  fixed  points  decrease  in  size  under  renormalization.  Such 
perturbations  are  said  to  be  “irrelevant.”  An  unstable  fixed  point  has  at  least 


46 


one  relevant  operator  which  grows  under  renormalization,  eventually  driving  the 
system  away  from  the  fixed  point  and  toward  a different  (stable)  fixed  point. 

The  set  of  initial  coupling  configurations  entering  H0  which  flow  to  a given 
fixed  point  defines  the  basin  of  attraction  or  domain  of  that  fixed  point.  The 
basin  of  attraction  of  a stable  fixed  point  has  a finite  measure  in  the  parameter 
space  of  the  model.  By  contrast  an  unstable  fixed  point  has  a domain  of  zero 
measure,  because  the  fixed  point  can  be  reached  only  if  the  parameters  are  tuned 
to  eliminate  all  relevant  perturbations. 

When  Vh  = 0,  the  HN’s  generated  by  R\  conserve  total  spin  5,  helicity 
h,  and  axial  charge  I (defined  in  Section  2.2).  All  the  fixed  points  identified 
below  obey  these  symmetries  and  possibly  involve  higher  symmetries  as  well.  We, 
therefore,  label  each  state  in  the  spectrum  of  any  of  these  fixed  points  according 
to  its  quantum  numbers  (5,  /,  h). 

3.1.2  Overview:  Particle-Hole-Symmetric  Phase  Diagram 

In  this  subsection  we  present  a preview  of  the  phase  diagram  of  the  particle- 
hole-symmetric  three-impurity  Kondo  model,  introducing  the  main  fixed  points. 
These  fixed  points  are  described  individually  in  Sections  3.2  and  3.3.  We  conclude 
the  chapter  with  a second,  more  detailed  description  of  the  phase  diagram. 

Figure  3.2  shows  a schematic  representation  of  a two-dimensional  projec- 
tion of  the  four-dimensional  parameter  space  of  the  particle-hole-symmetric  three- 
impurity  model.  The  horizontal  axis  of  Figure  3.2  measures  the  normalized  RKKY 
coupling 
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Figure  3.2:  Schematic  phase  diagram  for  the  three-impurity  Kondo  model  at 
particle-hole  symmetry.  Here  K is  a measure  of  the  RKKY  interaction,  and  a.  mea- 
sures the  relative  strengths  of  the  couplings  Joi  and  J\i  in  Eq.  (2.34).  Filled(hollow) 
circles  indicate  stable(unstable)  fixed  points.  The  heavy  broken  line  represents 
a line  of  unstable  fixed  points.  The  dotted  line  marks  the  border  between  the 
frustrated-Kondo  and  independent-impurity  domains 
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K = 


Jl o + J2  - A - 2 J,2 


'ii 


01 


‘■t  kky 


(3,5) 


+ 2,/?,  + 2,7=2  + Ii,?,  2 In  2(()„  J)2D  ’ 

The  vertical  axis  represents  the  relative  strength  of  the  exchange  couplings  which 
contribute  to  antiferromagnetic  RKKY  couplings,  as  parameterized  by 


a — 


Jm  — J 


12 


(3.6) 


Joi  + ^12 

With  these  definitions,  K ranges  from  -1/2  to  +1,  while  — 1 < a < +1.  Neither 
K nor  a depends  on  the  exchange  coupling  parameter  J.  This  allows  us  to  locate 
strong-coupling  fixed  points  on  a diagram  which  spans  only  the  finite- valued  initial 
coupling  configurations. 

The  phase  diagram  is  dominated  by  the  domains  of  three  stable  fixed  points. 
Throughout  the  region  where  the  RKKY  coupling  is  ferromagnetic  (K  > 0)  or 
weakly  antiferromagnetic  ( K < 0),  the  low-temperature  behavior  is  described  by  a 
strong-coupling  “independent-impurity”  fixed  point  at  which  the  three  impurities 
are  Kondo-screened  independently  of  each  other.  In  the  region  of  strong  antifer- 
romagnetic RKKY  couplings  ( K 0),  two  stable  non- Fermi-liquid  (NFL)  fixed 
points  exist:  the  “frustrated- Kondo”  for  positive  a and  the  “isospin  two-channel’ 
for  negative  a.  Unstable  fixed  points  mark  the  boundaries  between  these  stable 
fixed  points:  a single  unstable  NFL  fixed  point  separates  the  frustrated-Kondo 
and  independent-impurity  domains,  while  a line  of  unstable  fixed  points  surrounds 
the  stable  NFL  fixed  point  associated  with  negative  values  of  a.  Finally,  there  are 
two  unstable  Fermi-liquid  fixed  points  which  can  occur  at  the  strong  ferromagnetic 
RKKY  limit  and  a — 0 (but  distinguished  by  different  values  of  J0o  and  Jn). 
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3.2  Fermi-Liquid  Fixed  Points 

In  this  section  we  describe  the  fixed  points  which  represent  Fermi-liquid 
phases  of  our  three-impurity  model  and  identify  the  low-temperature  behavior 
associated  with  each  fixed  point. 

All  the  Fermi-liquid  fixed-point  spectra  map  onto  spectra  generated  by  a 
three-channel  free-electron  Hamiltonian 


(3J> 

where 

H^h=  £ + fUM,  (3-8) 

n=0  <r 

and  h serves  as  the  channel  index.  The  transformation  R\  generates  a different 
free-electron  fixed-point  spectrum  for  even  and  odd  N.  We  differentiate  the  Hamil- 
tonians which  generate  these  spectra  by  placing  a hat  over  the  even-N  Hamiltoni- 
ans. Thus  H' s generate  the  odd- A''  spectra  and  H's  generate  the  even -N  spectra 
of  the  free-electron  model. 

These  spectra  differ  due  to  the  different  symmetry  requirements  for  spectra 
generated  by  a tight-binding  Hamiltonian  acting  on  N + 1 sites.  When  N + 1 
is  odd  (i.e.,  N even),  there  are  necessarily  zero-energy  eigenvalues  in  addition  to 
both  positive  and  negative  energy  eigenvalues  of  H%.  On  the  other  hand  for  N + 1 
even  (N  odd),  there  are  no  zero-energy  eigenvalues  of  H°N.  This  argument  holds 
for  each  channel  of  the  free-electron  Hamiltonian. 
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Table  3.1:  The  Hamiltonian  and  symmetry  operators  of  the  free-electron  model 
for  odd  N . See  Table  3.3  for  limiting  values  of  the  coefficients  rjj  and  anj]  h = 
(3  — h ) mod  3. 


(JV+l)/2 


^N,h  ~ ^2  9]  ]\ho9jho  + hjhohjho) 

j=l  ° 

(3.9) 

(AT+l)/2 

Sjv,/i  = 5 ^2  U9]h^9^'  + hjhaaayhjha') 

j= 1 <r,cr' 

(3.10) 

(JV+l)/2 

•fv,/i  = 5 X]  ^2(9jhtr9jh<T  ~~  ^jhahjhc) 
j= 1 

(3.11) 

(JV+l)/2 

iw=  J2  - sJAm) 

j=l  /l 

(3.12) 

(tf+l)/2 

In  = E 'Y2^l)h[9fh]  ~ ^jh]9jhi) 

j= 1 h 

(3.13) 

(N+ 1)/2 

foha  = A"(yv_1)/4  a0j(^fca  + 

j=l 

(3.14) 

(JV+l)/2 

/lfca  = A-3(N-1)/4  E alj(Pjha  “ *»}fca) 

i=i 

(3.15) 
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Table  3.2:  Hamiltonian  and  symmetry  operators  of  the  free-electron  model  for  even 
N.  See  Table  3.3  for  limiting  values  of  the  coefficients  r)j  and  anj\  h = (3— h)  mod  3. 


^N,h  — X 

Sjv,h  = \ X 

cr, a' 

^N,h  = \ X 

a 

In  = X 

h 

In  = X 

h 

f0h,  = A-("-1)/4 


N/2 

Vodohirdjha  + Vj{9jh<j9jhcr  + ^jhcr^jha) 
3= 1 


N/2 

9oh<T^a,<r'90ha'  + ^2(9jha&<7,<r'  9jha'  + h^h(7^a,a' hjha') 

3= t 

N/2 

{9ohcr90hcr  ~ 1/2)  + ^^{9 j ha 9 j he  ~ ^'jha^'jhtr) 


3= 1 


N/2 


9oh]9ohi  + X 9jh]hfhl  ~ 9]hihjh]) 

3= 1 


N/2 


So/ijtfo/iT  + — 

j=l 

N/2 

Oi0090h<r  + X ( 9jh<r  + A jha) 


3=1 


N/2 


flha  = A 3(iV  1)/4  X 

3=1 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


(3.22) 
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Table  3.3:  Asymptotic  values  of  coefficients  entering  Tables  3.1  and  3.2  for  1 < 

j < N. 


r)j  = A(j'1} 
ctoj  = «oA^_1^2 


rjj  = A^  1 ^ with  r}0  = 0 
do,  = a0A^-1^2)/2 
= arA3^1/2)/2 


<*0  = 
«*i  = 


n 1/2 


1(1  -A"1) 

2(1  — A-3)l 1/2 
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Since  H°N  is  quadratic  it  can,  in  principle,  be  diagonalized  exactly.  Tables  3.1 
and  3.2  show  the  diagonal  representations  of  Hu  in  terms  of  particle  and  hole  op- 
erators, along  with  the  operators  which  generate  symmetries  of  the  fixed  point.  As 
N goes  to  infinity,  and  become  equivalent  to  the  fixed-point  Hamiltonians, 
H*  and  H*,  whose  asymptotic  spectral  properties  are  described  in  Table  3.3.  Note 
that  a spectrum  calculated  numerically  is  the  full  spectrum  corresponding  to  all 
possible  combinations  of  single-particle  excitations  gjha  and  consistent  with 
the  Pauli  exclusion  principle. 

3.2.1  Weak-Coupling  Fixed  Point 

The  weak-coupling  fixed  point  of  the  three-impurity  model  is  analogous  to 
the  weak-coupling  fixed  points  of  the  one-  and  two-impurity  Kondo  models.  It 
describes  the  behavior  of  the  three-impurity  model  when  the  conduction  band 
decouples  from  the  impurities.  Any  nonzero  exchange  coupling  drives  the  model 
away  from  weak  coupling  and  towards  one  of  the  other  fixed  points  of  the  model. 
For  example,  perturbations  in  J00  at  this  fixed  point  will  drive  the  model  to  the 
R — 0 fixed  point  (documented  below)  even  though  the  R = 0 fixed  point  is  itself 
unstable  to  perturbations  in  all  the  other  exchange  couplings. 

Except  for  an  additional  eight-fold  degeneracy  of  each  energy  level,  the  weak- 
coupling  fixed-point  spectra  map  the  free-electron  spectra  exactly, 

lim  ^H°Nh  = H*  (weak-coupling). 


(3.23) 
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The  additional  degeneracy  results  from  the  eight  possible  impurity  configurations 
at  this  fixed  point. 

As  with  all  the  other  Fermi-liquid  fixed-point  spectra  of  the  three-impurity 
Kondo  model,  the  addition  of  weak-to-moderate  potential  scattering  Vh  at  the 
weak-coupling  fixed  point  shifts  the  energy  of  each  state  by  an  amount  proportional 
to  its  charge,  Q = 2 Iz.  Near  the  Fermi-liquid  fixed  points,  potential  scattering 
represents  a marginal  operator  which  converts  each  particle-hole-symmetric  fixed 
point  into  a line  of  fixed  points. 

3.2.2  Independent-Impurity  Fixed  Point 

The  independent-impurity  fixed  point  is  analogous  to  the  strong-coupling 
fixed  point  of  the  one-impurity  Kondo  model  and  the  independent-impurity  strong- 
coupling fixed  point  of  the  two-impurity  model.  The  impurity  separation  is  effec- 
tively infinite,  as  the  Kondo  effect  screens  the  individual  impurities  separately. 
The  exchange  couplings  all  go  to  strong  coupling. 

The  independent-impurity  fixed-point  spectra  describe  three  noninteracting 
channels  of  conduction  electrons.  The  even-iteration  independent-impurity  spec- 
trum maps  the  odd-iteration  free-electron  spectrum, 

lim  ^H%h  = H*  (independent-impurity).  (3.24) 

°°  h 

The  ground  state  includes  no  impurity  or  electronic  degrees  of  freedom.  As  in 
the  spectra  calculated  by  Wilson  for  the  one-impurity  Kondo  model  discussed  in 
Section  1.2,  the  fact  that  the  even-iteration  spectra  generated  at  the  independent- 
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impurity  fixed  point  are  the  same  as  the  odd-iteration  spectra  generated  at  the 
weak-coupling  fixed  point  indicates  a 7r/2  phase  shift  in  the  electronic  wave  func- 
tions at  the  Fermi  surface  for  each  helicity  of  electron.  The  Kondo  effect  localizes 
one  electron  of  each  helicity  about  the  impurities,  thereby  “freezing  out  the  im- 
purity and  on-site  electronic  degrees  of  freedom. 

This  fixed  point  is  stable  to  perturbations  in  the  exchange  couplings.  As 
in  the  one-  and  two-impurity  Kondo  models,  the  exchange  couplings  of  each 
helicity-channel  tend  toward  strong  coupling  as  temperature  decreases  in  the  three- 
impurity  model  unless  prevented  by  other  processes. 

3.2.3  R = 0 Fixed  Point 

The  R = 0 fixed  point  describes  the  low-temperature  behavior  of  the  three- 
impurity  model  when  the  impurity  separation  equals  zero.  The  RKKY  coupling 
takes  its  maximum  ferromagnetic  value  causing  the  impurities  to  lock  into  a spin-| 
configuration.  Due  to  the  vanishing  of  the  structure  factors  gi(0)  and  #2(0)  [see 
Eq.  (2.5)],  only  the  helicity-0  electrons  couple  to  the  impurity  configuration.  The 
Joo  coupling  alone  goes  to  strong  coupling  and  a localized  helicity-0  conduction 
electron  partially  screens  the  impurity  configuration,  leaving  a residual  spin-1  mo- 
ment. As  a result,  the  fixed-point  spectra  incorporate  a 7t/2  phase  shift  in  the 
helicity-0  sector  of  the  excitation  spectra  and  no  phase  shift  in  the  h ^ 0 sector, 

lim  ( Y.  Kj.  + <o]  = (R  = 0)' 

N— >oo  " ’ 


(3.25) 
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This  fixed  point  is  unstable  to  perturbations  in  exchange  couplings  other 
than  J00,  which  drive  the  model  to  the  independent-impurity  fixed  point. 

3.2.4  Ju  (Strong-Coupling)  Fixed  Point 

This  fixed  point  does  not  reproduce  the  behavior  of  any  physically  realizable 
scenario  in  the  three-impurity  model.  At  the  J\\  fixed  point,  the  RKKY  interaction 
takes  its  maximum  ferromagnetic  value  just  as  in  the  R = 0 case,  but  here  only 
helicity-1  and  -2  electrons  couple  to  the  impurities. 

The  fixed-point  spectra  are  free-electron-like  with  a 7t/2  phase  shift  in  the 
nonzero- helicity  single-particle  excitations, 

lim  [ Y.  = »'  (Ju)-  (3  26) 

lit 2 

Due  to  the  strong  ferromagnetic  RKKY  coupling  a spin-|  configuration  forms  and 
is  partially  screened  to  result  in  a residual  spin-|  moment. 

This  fixed  point  is  unstable  with  respect  to  perturbations  in  Joo  and  Joi, 
which  drive  the  model  to  the  independent-impurity  fixed  point.  This  fixed  point 
is,  however,  stable  to  perturbations  in  Jn. 

3.3  Non-Fermi-Liquid  Fixed  Points 

The  non- Fermi- liquid  (NFL)  fixed  points  of  the  three-impurity  Kondo  Hamil- 
tonian all  occur  in  the  regime  of  moderate-to-strong,  antiferromagnetic  RKKY 
coupling.  Two  stable  NFL  fixed  points  occur  when  the  RKKY  coupling  takes  its 
maximum  antiferromagnetic  value,  while  unstable  NFL  fixed  points  occur  when 
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the  RKKY  coupling  and  the  Kondo  energy  D exp(l/ q0J)  — & characteristic  scale 
for  the  single  impurity  Kondo  effect  — become  of  the  same  order. 

3.3.1  Isospin  (Two-Channel)  Fixed  Point 

The  isospin  fixed  point  is  reached  whenever  the  initial  couplings  are  such  that 
the  RKKY  coupling  is  strongly  antiferromagnetic  and  J12  » Joi-  The  fixed  point 
spectra  for  even  and  odd  N are  listed  in  Table  3.4.  These  spectra  can  be  decom- 
posed into  products  of  free-electron  spectra  for  helicity-0  electrons  and  a spectrum 
which  cannot  be  described  in  terms  of  independent  free-electron-like  excitations. 
Such  a decomposition  can  be  explained  by  both  Joo  and  J01  renormalizing  to  zero 
at  the  fixed  point. 

The  helicity-0  spectrum  describes  a standard  Fermi  liquid.  More  interesting 
is  the  second,  non- Fermi-liquid  spectrum,  which  is  identical  for  even  and  odd 
iterations;  an  interesting  property  it  shares  with  the  single-impurity  two-channel 
Kondo  model.  Comparison  of  the  spectra  indicates  that  every  state  in  the  three- 
impurity  spectrum  has  a counterpart  of  equal  energy  in  the  two-channel  spectrum. 
The  quantum  numbers  of  equivalent  states  in  the  two  spectra  are  not  the  same. 
However,  one  finds  empirically  that  the  isospin  quantum  number  of  each  three- 
impurity  state  maps  onto  the  spin  of  the  two-channel  state.  In  particular,  the 
three-impurity  ground  state  is  an  isospin  doublet  — hence  our  name  for  this  fixed 
point.  In  fact,  as  we  show  in  Chapter  4,  one  can  exploit  a new  symmetry  present 
near  the  isospin  fixed  point  to  define  a complete  one-to-one  mapping  between  the 
isospin  spectrum  and  that  of  the  two-channel  fixed  point. 


58 


Table  3.4:  Low-lying  states  at  the  isospin  fixed  point.  Energies  E are  given  for  odd 
and  even  iterations.  The  energies  were  calculated  for  a discretization  parameter 
A = 3.0,  then  extrapolated  to  the  continuum  limit  by  multiplying  by  0.625,  a factor 
which  maps  the  first  excited  levels  of  the  A = 3.0  free-electron  spectra  onto  those 
for  A = 1.  States  are  labeled  by  their  spin  S , isospin  /,  and  helicity  A;  an  asterisk 
stands  for  one  state  each  of  helicity  1 and  2.  See  the  text  for  further  details. 


N even 

N odd 

E 

(5,  /,  h) 

E 

(S,  I,  h ) 

0 

(0,0,0)  (0,1,0) 

0 

(0,  J,0) 

55 

0.124 

(3,0,.) 

0.125 

(0,0,*)  (1,0,*) 

0.498 

(0,i,»)  (1,1,0) 

55 

(I  1 *) 

0.501 

(1,0,0) (i, 1,0) 

0.500 

(0,0,*)  (0,1,*) 

0.625 

(o,  J,  *)  (1,5,*) 

55 

(£>i°)  (ii*) 

0.628 

(5, 1,  *) 

55 

(1,0,0)  (1,1,0) 

0.997 

(0,|,0) 

55 

(§» 2>°) 

1.000 

(i,0,0)  (J,  1,0) 

0.629 

(0,1,*)  (1,1,*) 

55 

(1,0,.)  (1,1,.) 

55 

(ii*)  (§>§>*) 

55 

(1,0,0)  (§,1,0) 

1.001 


1.003 
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Table  3.5:  Low-lying  states  from  the  finite-size  spectrum  at  the  frustrated-Kondo 
fixed  point.  See  Table  3.4  for  an  explanation  of  the  entries. 


N even 

E 

(S, 

I,h) 

0 

(0,0,0) 

x2 

0.100 

(I  1 *) 

0.200 

(0,1,*) 

(1,0,*) 

0.299 

(1,1,  o) 

x2 

0.499 

(0,0,*) 

x2 

0.599 

a i *) 

0.601 

(i.i.*) 

0.604 

(1,1, o) 

x2 

0.699 

(0,1,0) 

(1,0,0) 

0.700 

(0,1,0) 

(1,0,0) 

0.701 

(ill*) 

(i.i*) 

0.704 

(0,1,*) 

(1,0,*) 

0.804 

N odd 

E (5,  /,  h) 


0 

(oq,*) 

(i 

0,*) 

0.400 

d.i.*) 

(1, 

i*) 

0.500 

(0,1,0) 

x2 

(io,o) 

x2 

0.501 

(0,1.*) 

(i 

0,*) 

0.603 

(1,1,0) 

x2 

?) 

(1,1,0) 

x2 

0.799 

(0,f,0) 

(^ 

V2’ 

0,0) 

0.800 

(0,|,0) 

\2' 

0,0) 

0.904 

(1,1,0) 

(li 

i.o) 

0.905 

(1,1,0) 

(1, 

2’0) 

0.913 

(i  1)  *) 

(1, 

i*) 

1.00 

This  fixed  point  is  stable  with  respect  to  perturbations  in  the  exchange  cou- 
plings. Particle-hole  asymmetry  drives  the  model  to  a Fermi- liquid  fixed  point,  a 
point  we  will  revisit  in  Chapter  4. 

3.3.2  Frustrated-Kondo  Fixed  Point 

We  present  the  low-energy  spectra  of  the  frustrated-Kondo  fixed  point  in 
Table  3.5.  We  cannot  map  the  frustrated-Kondo  spectra  onto  any  known  spectra 
and  believe  they  correspond  to  a completely  novel  NFL  behavior. 
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The  characteristics  which  identify  these  spectra  as  non-Fermi  liquid  are  as 
follows:  (i)  There  are  two  (0, 0, 0)  ground  states  in  the  iV-even  spectrum,  (ii)  States 
which  in  a Fermi-liquid  spectrum  would  incorporate  more  than  one  electron  or  hole 
excitation  have  a lower  energy  than  states  which  would  incorporate  a single  electron 
or  hole  excitation  (based  on  their  charge  quantum  number),  (iii)  A gap  exists  in 
the  N-odd  spectrum  at  E — 0.400.  In  addition,  an  analysis  of  the  approach  of 
the  energy  levels  of  H n to  their  fixed-point  values  with  increasing  N allows  one 
to  extract  the  exponent  A entering  the  Tx  temperature  variation  of  the  leading 
irrelevant  operator (s).  The  observed  value  A = 0.21±0.03  is  very  different  from 
the  value  A = 1.0  characteristic  of  a Fermi  liquid. 

This  fixed  point  is  stable  with  respect  to  perturbations  in  all  exchange  in- 
teractions and  is  not  immediately  destroyed  by  particle-hole  asymmetry.  Instead, 
weak-to-moderate  potential  scattering  shifts  the  energy  levels  of  each  state  by  an 
amount  proportional  to  the  charge  of  the  state,  A Ey  oc  VIZ.  However,  strong 
potential  scattering  drives  the  model  to  a weak-coupling  fixed  point.  A regime  of 
transitional  behavior  occurs  for  intermediate  values  of  14. 

3.4  Particle-Hole-Symmetric  Phase  Diagram 

In  this  section  we  present  a more  detailed  discussion  of  the  phase  diagram  of 
the  three-impurity  model.  We  start  by  explaining  the  choice  of  the  parameters  K 
and  a as  axes  for  Figure  3.2.  We  then  present  a detailed  discussion  of  the  interplay 
of  the  Kondo  effect  and  RKKY  interaction  in  the  phase  diagram. 

The  quantities  K and  a as  defined  in  Eqs.  (3.5)  and  (3.6)  are  natural  choices 
for  construction  of  the  phase  diagram,  because  they  describe  the  competitions 
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which  dominate  the  behavior  of  the  three-impurity  Kondo  model.  K parameter- 
izes the  competition  between  magnetic  ordering  of  the  impurities  and  quenching 
through  the  Kondo  effect,  while  a measures  the  relative  strength  of  the  interac- 
tions which  compete  to  drive  the  system  to  one  or  the  other  of  the  stable  non- 
Fermi-liquid  fixed  points.  We  find  the  value  of  K correlates  strongly  with  the 
low-temperature  impurity-spin  configuration  as  measured  by  the  spin  correlation 
function  ± • S j),  which  varies  smoothly  from  1/4  through  0 to  -1/4  as  K 

goes  from  +1  through  0 to  -1/2,  regardless  of  the  fixed  point  the  system  eventually 
flows  to. 

The  parameter  a is  a good  discriminator  in  the  antiferromagnetic  RKKY 
regime.  The  model  will  flow  to  either  the  frustrated-Kondo  or  the  independent- 
impurity  fixed  point  when  a > 0,  and  to  either  the  isospin  or  independent-impurity 
fixed  point  when  a < 0. 

The  variation  of  behavior  along  the  Y-axis  of  the  phase  diagram  shows  some 
similarities  with  the  two-impurity  model.  The  domain  of  the  independent-impurity 
fixed  point  encompasses  a broad  range  of  parameter  space  in  which  the  RKKY 
coupling  is  strongly  ferromagnetic  to  weakly  antiferromagnetic.  The  R = 0 and 
Jn  fixed  point  represent  the  only  exceptions.  As  we  have  already  described,  these 
fixed  points  are  reproduced  only  when  the  RKKY  coupling  takes  its  maximum 
ferromagnetic  value.  They  have  a direct  counterpart  in  the  underscreened  5=1 
fixed  point  of  the  two-impurity  model  [31]. 

The  two-  and  three-impurity  models  behave  very  differently  in  the  limit  of 
strong  antiferromagnetic  RKKY  couplings.  In  the  former  model  the  impurities 
can  lock  into  a spin-singlet,  causing  the  electrons  to  experience  pure  potential 
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scattering  from  the  impurities.  In  the  three-impurity  case,  the  impurities  form  a 
frustrated  spin-doublet,  and  one  of  two  novel  regimes  is  entered.  This  confirms 
our  hypothesis  that  frustration  provides  an  alternative  route  to  non-Fermi-liquid 
behavior. 

Unstable  NFL  fixed  points  lie  on  the  boundaries  between  the  NFL  domains 
and  the  independent-impurity  domain.  Non-Fermi  liquid  behavior  occurs  on  these 
boundaries  due  to  the  noncommutivity  of  the  Kondo  effect  and  RKKY  interaction 
in  the  regime  of  antiferromagnetic  RKKY. 

In  the  two-impurity  model,  there  is  also  an  unstable  NFL  fixed  point  sepa- 
rating the  stable  regimes  of  the  model.  An  important  difference,  though,  is  that 
the  position  of  this  fixed  point  is  universal  in  the  two-impurity  case  (it  occurs  at 
Irkky/TK  = -2.4  for  all  values  of  J and  R).  By  contrast,  there  is  no  simple  scal- 
ing of  Irkky  (or  K)  by  TK  which  eliminates  all  J dependence  from  the  size  of  the 
three  stable  domains  of  the  three-impurity  model.  Since  the  Kondo  temperature 
increases  as  exp(l/p0^)  while  the  RKKY  interaction  is  quadratic  in  J,  increasing 
| J\  will  increase  the  size  of  the  independent-impurity  domain  at  the  expense  of  the 
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CHAPTER  4 
ANALYSIS  OF  RESULTS 


This  chapter  presents  further  analysis  of  the  stable  fixed  points  of  our  three- 
impurity  Kondo  model.  The  first  section  describes  a flow  analysis  leading  to  the 
determination  of  the  Wilson  ratio  for  various  regions  within  the  domain  of  the 
independent-impurity  fixed  point.  The  second  section  presents  an  analysis  of  the 
isospin  fixed  point.  We  describe  a minimal  model  which  reproduces  the  fixed 
point  spectrum.  We  also  derive  the  low-temperature  properties  using  a mapping 
between  the  isospin  fixed  point  and  the  non- Fermi-liquid  fixed  point  of  the  two- 
channel  spin-  \ Kondo  model.  The  final  section  deals  with  the  frustrated-Kondo 
fixed  point,  where  we  identify  two  symmetries  which  are  present  at  the  fixed  point. 

4.1  Independent-Impurity  Domain 

In  this  section  we  analyze  the  properties  of  the  three-impurity  model  in  the 
basin  of  attraction  of  the  independent-impurity  fixed  point  using  a perturbative 
expansion  about  the  fixed  point  in  terms  of  an  effective  Hamiltonian  whose  coupling 
parameters  may  be  determined  from  the  low-energy  spectrum  of  near  the  fixed 
point. 

One  can  calculate  thermodynamic  properties  from  the  spectra  generated  by 
Hn  in  two  ways:  either  a “brute-force”  calculation  performed  by  summing  the  con- 
tribution from  each  eigenstate  retained  in  the  calculation  (as  discussed  by  Oliveira 
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and  Oliveira  [50])  or  a perturbative  formalism  developed  by  Wilson  [72].  The 
accuracy  of  the  brute-force  approach  is  determined  by  the  number  of  eigenstates 
retained  from  one  iteration  to  the  next.  Suppose  that  one  is  able  to  retain  all 
eigenstates  with  scaled  energies  (measured  with  respect  to  the  ground  state)  up  to 
Ec,  where  the  smallest  splitting  is  by  definition  of  order  unity.  Since  the  character- 
istic temperature  at  iteration  N is  [see  Eq.  (4.9)]  X/v  ~ A ~N^2D/kB,  after  scaling 
the  thermal  scale  is  of  order  unity.  As  a result  the  finite  value  of  Ec  introduces 
errors  into  the  brute-force  method  of  evaluating  thermodynamic  sums  at  the  level 
of  e~Ec. 

In  the  one-impurity  Kondo  model  it  is  possible  to  work  with  Ec  > 20.  How- 
ever, the  larger  Hilbert  space  of  the  three-impurity  model  prevents  calculation  of 
eigenstates  of  HN  with  Ec  much  greater  than  5 and  effectively  rules  out  brute- force 
calculation  of  thermodynamic  properties  for  this  model. 

Near  a Fermi-liquid  fixed  point,  such  as  the  independent-impurity  fixed  point, 
the  problem  of  working  with  a basis  restricted  by  truncation  can  be  circumvented 
through  the  use  of  an  effective  Hamiltonian  which  generates  the  entire  fixed-point 
spectrum.  This  allows  one  to  calculate  thermodynamic  properties  perturbatively  in 
couplings  of  operators  describing  deviations  from  the  fixed-point  Hamiltonian  [72]. 
The  numerical  values  of  these  couplings  (or  at  least  those  couplings  which  make  the 
greatest  contribution  to  the  thermodynamics)  can  be  extracted  from  the  details  of 
the  approach  of  each  eigenvalue  of  Hn  to  its  fixed-point  value. 

The  effective  Hamiltonian  consists  of  a finit e-N  approximation  to  the  fixed- 
point  (i.e. , N — > oo)  Hamiltonian  plus  a term  which  describes  deviations  from  the 
fixed  point, 
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HeNff  = H*n  + AHn.  (4.1) 

Both  H*n  and  A HN  act  on  a chain  of  N+l  sites.  H*N  has  the  symmetries  of  the 
fixed  point  Hamiltonian,  H*  — limAr_oo  H*N.  The  deviation  A contains  all 
allowed  perturbations  about  the  fixed  point.  It  can  be  written  in  the  form 

AHN  = YJ^N-l)/28Hi.  (4.2) 

i 

Here  SHi  is  an  operator  composed  from  any  of  the  degrees  of  freedom  which  remain 
active  at  the  fixed  point.  At  the  independent-impurity  fixed  point,  this  means  any 
fnhc  with  n > 0.  (The  /0/l<7  and  impurity  degrees  of  freedom  are  “frozen  out”  in 
the  strong-coupling  limit.) 

Each  perturbation  8Hi  is  constrained  to  obey  all  the  symmetries  of  the  initial 
three-impurity  problem,  but  may  break  the  higher  symmetries  of  the  fixed  point 
about  which  one  is  expanding.  The  factor  A^-1^2  which  enters  Eq.  (4.2)  simply 
reflects  the  rescaling  which  takes  place  at  each  iteration  of  the  recursion  relation 
for  Hn  [Eq.  (3.1)].  The  coefficient  u>i  represents  the  weight  factor  for  5Ht  in  A H^. 
The  value  of  u \ depends  on  the  bare  couplings  J*  and  V*  entering  Ho- 

Very  close  to  a fixed  point,  the  deviation  A Hn  must  by  definition  be  very 
small.  Thus  it  is  possible  to  expand  the  renormalization-group  transformation 
equation  [Eq.  (3.1)]  about  H*N  to  obtain  a linearized  transformation  C*A: 

A Hn+2  = Rl[H*N  + AHN]-Rl[H*N\ 


£(H'n)&hn 


(4.3) 
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with  C(h*n)  = C(H*)  = for  large  N-  Since  C*N  is  a linear  operator,  it  generates 
eigenoperators  0*  with  eigenvalues  A*.  Through  repeated  application  of  C\,  one 
sees  AHn  can  be  written 


(4.4) 

i 

Here  the  fact  approximates  R\  rather  than  R\  leads  to  the  factor  of  | in 

the  power  of  A*.  The  operators  0*  are  formed  not  from  the  original  degrees  of 
freedom  fnh(T,  but  from  the  particle  and  hole  operators  which  diagonalize  H*N.  In 
the  domain  of  the  independent-impurity  fixed  point,  these  are  the  gjha  s and  hjha  s 
entering 

(N- 1)/2 

hn=  L VjWjh^jhc  + h)h(Thjha]  (independent-impurity).  (4.5) 
j= o 

We  define  t)*  and  the  particle  and  hole  operators  in  our  discussion  of  the  free- 
electron  Hamiltonian  in  Section  3.2  and  Tables  3.1  and  3.2. 

It  is  important  to  note  that,  at  least  for  large  N,  O*  is  independent  of  N. 
Then  Eq.  (4.4)  implies  0*  is  relevant  (grows  in  importance  with  increasing  N) 
if  the  associated  eigenvalue  A*  is  greater  than  one  and  is  irrelevant  (diminishes 
in  importance)  if  the  associated  eigenvalue  is  less  than  one.  In  the  vicinity  of  a 
stable  fixed  point,  all  the  operators  must  be  irrelevant.  Moreover,  for  large  N one 
need  only  consider  the  effect  of  the  leading  or  least  irrelevant  operators,  i.e.,  those 
having  the  largest  (common)  value  of  A*. 

In  practice  one  constructs  the  perturbations  5 Hi  from  the  operators  fnha- 
Each  5Hi  can  then  be  expanded  in  terms  of  the  particle-hole  operators  which 


67 


diagonalize  H*N  by  expanding  the  fnhc' s.  (Tables  3.1  and  3.2  present  expansions  for 
f0h<j  and  fih(r  in  terms  of  the  particle  and  hole  operators  which  diagonalize  the  free- 
electron  fixed  point.)  One  can  then  extract  the  eigenvalue  of  each  eigenoperator 
from  the  relation  implied  by  Eqs.  (4.2)  and  (4.4),  namely, 


At  the  independent-impurity  fixed  point,  the  expansion  of  f\ha  incorporates  a 
factor  of  A-l^-1)/4  and  that  of  f2ha  incorporates  a factor  of  A-3(iV-1)/4.  The  fnha' s 
with  n > 2 incorporate  even  smaller  prefactors.  Therefore,  the  leading  irrelevant 
operators  (LIOs)  generated  near  the  independent-impurity  fixed  point  are  those 
formed  from  the  fewest  fih(r's  and  f2hc  s.  Operators  which  involve  a single  fnha  are 
ruled  out  by  charge  conservation,  while  those  involving  flhafnh'a1  violate  particle- 
hole  symmetry.  Thus  the  LIOs  generated  near  the  independent-impurity  fixed 
point,  all  having  A*  = A-1,  involve  four  ,/W’s  or  an  fiha  and  an  f2ha-  We  present 
a set  of  interactions  which  generate  a complete  set  of  LIOs  valid  in  the  basin  of 
attraction  of  the  independent-impurity  fixed  point  in  Table  4.1.  Here  the  effective 
Hamiltonian  which  approximates  H ^ takes  the  form 


(4.6) 


HeJf  = h*n  + ^2  CiO* , 


(4.7) 


with  Ci  = A-W-Md. 


We  label  the  LIOs  in  Table  4.1  based  on  the  type  of  interaction  which  gener- 
ates them.  A form  of  Coulomb  interaction  between  the  site-1  electrons  generates  all 
the  operators  labeled  17,,  an  exchange  interaction  involving  site-1  electrons  gener- 
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Table  4.1:  Definitions  of  the  leading  irrelevant  operators  for  the  independent- 
impurity  fixed  point.  The  operators  n/,  = and  — fihafiha  ■ 


0Uo  = A(iV  1}(-n0  + 2n0Tn0|  + 1) 

= A^(no-l)2 

Oui  = A(iV_1^(— 5ni  — 5ri2  + 2ni|ni|  + 2n2|n2|  + 2ni|?7,2t  + 2nijri2| 

+ 6n1Tn2i  + 6n2tnu  + + 4/i2T  /lil/m/iiT  + 8) 

0Uoi  = A(iV_1)[— 2n0  - ni  - n2  + nQ{nx  + n2)  + Zfloifloifmfni 

+2/io|/io|/lll/l2T  + 2/ii|/i2|/lO|/lOT  + 2/i2|/ll|/l0i/l0T  + 2] 

Oj01  = A(Af_1)(n0tni|  - n0|ni|  + n0|n2i  - not^i  + «o|^u  ~ wiTnox 

+ no|n2j  — n2Tn0j.  — 2/Ioi/^i/ioi/ht  — 2/10|/i2|/ioi/i2T 
— 2/ii|/ioj,/u/iot  ~ 2/i2|/ioj./i2i/ior) 

Ou'u  = A(Af_1)(-ni  -n2  + ^n^rin  + 4n2|n2|  - 2n1Tn2|  - 2 nun2j 

+ 2/i+iT/i2l/lli/l2T  + 2/i2|/lli/l2|/llT  + 4) 

Ot/Jm  = A(jV_1)(//oT/i2|/ll|/llT  + /10T/111/121/12T  + /l2T/l0|/ll|/llT 

T/iiT^ioi/m/m  + fn^fnifoifn]  + /iit/iii/i°l/i2T  + fn]f\2i  fulfil 

Oth  = + for  a = 0,1, 2 
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ates  all  operators  labeled  by  an  Jj,  and  an  interaction  where  electrons  hop  between 
site  1 and  2 generates  the  operators  labeled  U. 

We  calculate  the  impurity  contribution  to  the  susceptibility  of  the  three- 
impurity  model  in  the  independent-impurity  domain  by  substituting  Hj^  for  H n 
in  the  definition  of  the  impurity  susceptibility: 


kBTXimP(T) 


TrSzN2exy(-PHN)  TtSzn2  exp(-pH°N) 


(4.8) 


Trexp (-/3Hn)  TV  exp 

Here  the  second  term  subtracts  the  contribution  of  the  noninteracting  conduc- 
tion band  in  the  absence  of  the  impurity  configuration.  We  give  the  free-electron 
Hamiltonian  H%  and  the  free-electron  spin  operators  in  Tables  3.1  and  3.2, 
while 


/?=5£?(1  + A-‘)A^V-1,/2.  (4.9) 

Evaluation  of  the  traces  in  Eq.  (4.8)  is  simplified  by  using  the  values  of  various 
sums  calculated  in  Section  V of  Krishna- murthy  et  al  [36]. 

After  performing  the  traces,  perturbatively,  to  first-order  in  the  5 Hi,  using 
the  basis  of  particle  and  hole  excitations  which  diagonalize  HN  and  HQN,  we  find 


Ximp 


(g^s)2  2 

D 1 + A-1 


( Uq  + 6f/i  — 4Jq  + 6C/n) 


— (to  + + ^2) 


a0ai 

In  A J ’ 


(4.10) 


where  ao  and  a.\  are  defined  in  Table  3.3. 
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Similarly,  we  derive  the  impurity  contribution  to  the  Helmholtz  free  energy 
for  the  model  in  the  independent-impurity  domain  by  substituting  HeNff  for  HN  in 
the  expansion  of  the  free  energy: 

Fimv(T)  = -kBT  {\n[Tr  eM~PHN)}  - lnffrexp  (-0H°N)]}  . (4.11) 


Here  we  find 


FimP{T)  = ^ + t\+  t2) 


2q0Qi 

.I+A-IJ  3 In  A 


-7T 


(4.12) 


In  order  to  compute  the  numerical  values  of  thermodynamic  properties,  one 
needs  to  know  the  coefficients  Cj.  These  can  be  obtained  by  equating  energy 
splittings  of  low-lying  states  calculated  using  first-order  perturbation  theory  in  the 
8 Hi's  to  the  same  splittings  in  the  spectrum  generated  by  HN, 


(m\HeJf\m)  — ( m'\HeNff\m ')  - £W,m  — -EW.m' • (4-13) 

Here  E^<rn  and  Es,m'  are  the  numerical  eigenvalues  at  iteration  N of  the  states 
which  evolve  into  the  states  \m)  and  \m')  as  N — > oo.  One  solves  for  the  Cj’s  by 
inverting  a linearly  independent  set  of  equations  of  the  form  of  Eq.  (4.13). 

Based  on  the  expansions  for  Ximp  given  in  Eq.  (4.10)  and  the  expansion  for 
the  specific  heat  derived  from  the  expansion  in  Eq.  (4.12), 

c,mp  = , (4.14) 

we  calculate  an  expression  for  the  Wilson  ratio: 
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Figure  4.1:  Wilson  ratio  Rw  versus  scaled  RKKY  coupling  K [as  defined  in 
Eq.  (3.5)]  for  the  three-impurity  Kondo  model  defined  in  Eqs.  (2.28)-(2.32)  and 
(2.34),  with  A = 3,  V0  = Vi  = 0,  J0 o = Ju  = Ju,  and  various  values  of  J0 1.  Here 
the  dashed  line  represents  the  location  of  the  boundary  between  the  independent- 
impurity  and  frustrated- Kondo  domains. 


47t2  k%  Ximp  _ 1 _ {Up  + Ui  ~ 4Jqi  - Un) 

3 {gHB)2Cirnp  (to  + 2ti)  ailnA 


(4.15) 


Here  we  use  the  fact  that  t\  = t2  to  simplify  the  expression. 

We  present  a plot  of  the  Wilson  ratio  versus  the  RKKY  coupling  in  Fig.  4.1. 
This  set  of  data  was  obtained  by  setting  Jqo  = J\  1 = Ju  and  varying  Jqi-  These 
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results  and  others  calculated  throughout  the  domain  of  the  independent-impurity 
fixed  point  confirm  our  interpretation  of  the  three  Fermi-liquid  fixed  points  iden- 
tified in  Chapter  3.  For  runs  representing  infinite  impurity  separation  (all  Jt’s 
equal)  the  Wilson  ratio  equals  2,  as  one  would  expect  for  independently  quenched 
impurities.  Near  the  R = 0 fixed  point  the  Wilson  ratio  goes  to  2.27  ± 0.20  which 
approaches  the  value  2.0  calculated  using  conformal  field  theory  for  the  Wilson 
ratio  for  a single-channel  spin-|  Kondo  model  [2].  Near  the  J\\  fixed  point  the 
Wilson  ratio  goes  to  2.56  ± 0.20  which  approximates  8/3  derived  from  confomal 
field  theory  for  a two-channel  spin-|  Kondo  model  [2].  In  each  case,  the  value  of 
the  Wilson  ratio  confirms  our  picture  of  the  low-temperature  physics  at  the  fixed 
point.  These  results  also  validate  the  accuracy  of  our  numerical  calculations. 

4.2  Isospin  (Two-Channel)  Domain 

In  this  section  we  present  our  analysis  of  the  isospin  fixed  point  and  the 
properties  of  the  three-impurity  model  in  the  basin  of  attraction  of  the  isospin 
fixed  point.  We  formulate  a minimal  (simplified)  model  which  approximates  the 
full  model  near  the  isospin  fixed  point  and  reproduces  the  isospin  spectra  under 
renormalization.  We  argue  that  both  the  weak-  and  strong-coupling  fixed  points  of 
this  model  are  unstable  and  thus  the  isospin  fixed  point  occurs  for  finite  exchange. 
Finally  we  formulate  a mapping  between  the  spectra  of  the  isospin  fixed  point  and 
the  non-Fermi-liquid  fixed  point  of  the  two-channel  S = 1/2  Kondo  model. 
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4.2.1  Minimal  Model 

We  formulate  a minimal  model  which  describes  the  three-impurity  Kondo 
model  in  the  vicinity  of  the  isospin  fixed  point  by  dropping  the  helicity-0  conduc- 
tion electrons  from  the  problem,  since  these  electrons  decouple  from  the  impurity 
configuration  at  the  fixed  point.  We  also  restrict  the  impurities  to  spin-|  configura- 
tions, since  the  isospin  fixed  point  is  only  reproduced  for  strong-antiferromagnetic 
RKKY  couplings.  With  these  simplifications,  the  Hamiltonian  takes  the  form 


Hm/D  = HJD  + Jn So  • (sn  + s22)  + Ji2(S!  • s21  + S2  • s12).  (4.16) 

Here  the  impurity  spin  operators  act  only  on  the  spin-|  configurations  and  we 
define  s hh>  in  Eq.  (2.35).  The  conduction  band  represented  by  Hc  includes  only 
nonzero-helicity  electrons. 

We  find,  using  the  numerical  renormalization-group  technique  described  in 
Chapter  2,  that  for  all  Jn  ± 0 and  Jn  > 0 the  stable  fixed-point  spectrum  of  the 
minimal  model  reproduces  the  non- Fermi-liquid  part  of  the  spectrum  generated 
by  the  full  model  at  the  isospin  fixed  point.  This  confirms  our  conjecture  that 
the  impurity  configuration  is  locked  in  a frustrated  spin-doublet  and  the  helicity-0 
electrons  decouple  from  the  impurity  configuration  at  the  isospin  fixed  point. 

We  can  gain  insights  into  the  minimal  model  by  considering  its  atomic  limit, 
i.e.,  the  Hamiltonian  H0  constructed  from  the  interaction  term  in  Eq.  (4.16).  For 
all  —0.83  < Jn/\Ju\  < 4.8,  the  ground  states  of  Hq  are 
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5 (/o'iTl2  i>  - /ou|2  t>  + At  I1  i)  - /oVl1  T>)  V'  = -l/2)  (4.17) 


These  ground  states  have  the  same  quantum  numbers  as  the  ground  states  of 
the  isospin  spectrum  when  the  contribution  of  the  helicity-0  electrons  has  been 
removed,  namely,  S = 0,  / = |,/z  = ±|,h  = 0.  These  states  act  as  an  effective 
residual  “isospin  doublet.” 

Based  on  an  analogy  with  the  argument  of  Nozieres  and  Blandin  [49]  con- 
cerning the  stability  of  the  strong-coupling  fixed  point  of  the  two-channel  Kondo 
model,  we  conclude  that  for  large  values  of  Jn  and  J\2  the  ground  states  of  Ho  will 
be  overscreened  in  the  isospin  channel  by  electrons  at  site  1,  and  that  this  process 
should  continue  indefinitely  for  sites  further  along  the  chain.  Although  the  analogy 
is  not  perfect,  because  there  are  two  independent  couplings  in  the  present  model, 
this  argument  suggests  that  both  J\\  and  J\i  will  flow  to  finite  values,  thereby 
explaining  the  non- Fermi- liquid  behavior  of  this  fixed  point. 


and 


\ (/o1t/o2t/ohI2  i)  + /oit/o1|/o2iI2  T)  + /0IT/02T/021II  D 

(Iz  = +1/2). 


(4.18) 
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4.2.2  Mapping  to  Two-Channel  Fixed  Point 

The  relationship  between  the  isospin  fixed  point  and  the  non- Fermi- liquid 
fixed  point  of  the  two-channel  Kondo  model  can  be  made  explicit  by  exploiting 
an  additional  symmetry  present  in  the  minimal  model  but  not  in  the  full  three- 
impurity  Hamiltonian.  This  symmetry  leads  to  a new  conserved  quantum  number, 


£‘  = i£ 


£(/£/„!,,  - fLU 2a)  + IM(M  - |2<7)(2<7| 


_n=0 


(4.19) 


Here  Lz  can  be  combined  with  the  ladder  operators  L+  = X!a(X^=o  fnicfn2a  + 
|lcr)(2cr|)  and  L~  = (L+)f  to  form  a vector  operator  L,  whose  components  obey 
angular  momentum  commutation  relations.  Conservation  of  Lz  implies,  but  is  a 
stronger  condition  than,  conservation  of  helicity,  where  the  helicity  of  a state  is 
h — Lz  mod  3. 

Use  of  Lz  allows  us  to  establish  a one-to-one  mapping  between  the  spectra 
of  the  isospin  and  two-channel  fixed  points.  Specifically,  the  isospin,  spin,  and  L 
degrees  of  freedom  at  the  three-impurity  fixed  point  map  onto  the  spin,  flavor,  and 
isospin  degrees  of  freedom  (respectively)  at  the  two-channel  fixed  point: 


three-impurity  < 


I 

S 

L 


S 

J/ 

I 


> two-channel. 


(4.20) 


This  mapping  can  be  accomplished  at  the  level  of  the  creation  and  annihilation 
operators  via  the  transformations 
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fnl]  * * 9n\ T 

/nil  + * Qnl]  (4'21) 

fn  2T  * * (-I)n+15n21 

fn2l  * * (-1)^11- 

Here  gnca  annihilates  an  electron  at  site  n in  channel  c with  spin  z-component  a 

in  the  tight-binding  representation  of  the  two-channel  Kondo  model.  The  factor 
of  (— l)n+1  ensures  the  particle-hole  symmetry  of  both  fixed  points. 

The  mappings  in  Eqs.  (4.20)  and  (4.21)  allow  us  to  deduce  the  physical 
properties  of  the  three-impurity  Kondo  model  near  its  isospin  fixed  point  from  the 
properties  of  the  two-channel  Kondo  model.  In  particular,  the  impurity  contribu- 
tion to  the  specific  heat  will  be  [19] 

Cimp  °c  — XTnT.  (4.22) 

The  spin  susceptibility  at  the  isospin  fixed  point  will  have  the  same  temperature 
dependence  as  the  flavor  susceptibility  Ximpj  of  the  two-channel  model  at  its  fixed 
point. 

The  mappings  also  allows  us  to  understand  the  instability  of  the  isospin  fixed 
point  with  respect  to  particle-hole  asymmetry.  Such  an  asymmetry  maps  onto  a 
spin  asymmetry  in  the  two-channel  Kondo  model,  i.e.,  the  effect  of  a magnetic 
field.  It  has  long  been  known  that  a magnetic  field  destroys  the  non-Fermi-liquid 
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behavior  of  the  two-channel  Kondo  model  and  causes  a crossover  to  Fermi-liquid 
low-temperature  properties  [49]. 

4.3  Frustrated  Kondo  Domain 

This  section  documents  our  analysis  of  the  frustrated- Kondo  fixed  point.  As 
noted  in  Chapter  3,  the  frustrated- Kondo  fixed-point  spectra  do  not  match  any 
spectra  previously  known  from  other  quantum  impurity  models.  Thus  there  exists 
no  simple  mapping  which  allows  us  to  determine  the  properties  at  the  fixed  point. 
As  mentioned  in  Section  4.2  the  constraints  of  computer  memory  and  time  prevent 
brute-force  calculations  of  thermodynamic  properties  of  the  three-impurity  Kondo 
model.  Also,  the  perturbative  technique  we  employed  to  calculate  the  properties  of 
the  independent-impurity  fixed  point  cannot  be  applied  to  non-Fermi-liquid  fixed 
points.  We  therefore  limit  our  analysis  of  the  frustrated- Kondo  fixed  point  to  the 
discussion  of  a minimal  model  which  reproduces  the  fixed-point  spectra. 

Our  minimal  model  for  the  frustrated- Kondo  fixed  point  restricts  the  impu- 
rities to  frustrated  spin-4  configurations  and  sets  J12  = 0: 

Hfk/D  — HJ D + J00S0  • Soo  + <AiS  • (sn  + S22)  (4-23) 

+J01S1  • (S10  + S02)  + J01S2  • (soi  + S2o)- 

Here  we  use  the  same  impurity  and  conduction  electron  spin  operators  as  in  our 
minimal  model  for  the  isospin  fixed  point.  This  model  reproduces  the  low-energy 
spectra  of  the  frustrated- Kondo  fixed  point  in  the  limit  N — ♦ 00. 
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This  model  has  a higher  symmetry  than  the  full  three-impurity  model,  which 
leads  to  a new  conserved  quantum  number 


m*  = y. 


E(/nl*/»l«r  - fl2Jn2c)  + £ &M(M 


71=0 


(4.24) 


Note  that  aside  from  a factor  of  two,  Mz  differs  from  the  quantity  L~  only  in 
the  treatment  of  the  impurity  configuration.  Again,  the  new  quantum  number 
supersedes  helicity  so  that  we  can  label  each  eigenstate  of  the  minimal  model  by 
its  total  spin,  axial  charge,  and  Mz  value.  Note  that  Mz  always  takes  an  odd 
integer  value. 

The  spectra  generated  at  the  frustrated-Kondo  fixed  point  exhibit  two  addi- 
tional symmetries: 

1)  Each  states  with  quantum  numbers  (S,  /,  Mz)  where  S ^ I is  degenerate 
with  a state  in  which  the  spin  and  axial  charge  quantum  numbers  are  swapped, 
i.e,  E(S,  I,  Mz)  = E(I,  S,  Mz). 

2)  For  every  state  (S, 1,  Mz ) there  is  a degenerate  state  with  quantum  num- 
bers (5, 7,  —Mz). 

At  present  we  do  not  know  the  significance  of  these  degeneracies. 


CHAPTER  5 
CONCLUSION 


We  conclude  this  dissertation  by  comparing  our  initial  expectations  for  the 
three-impurity  Kondo  model  with  the  low-temperature  behaviors  we  have  actually 
found.  We  discuss  some  of  the  implications  of  these  results  and  suggest  some 
avenues  for  future  research. 


5.1  Summary 

As  we  stated  in  Section  1.5  we  expected  the  behavior  of  our  three- impurity 
Kondo  model  to  be  similar  to  those  seen  in  the  two-impurity  model  for  regions 
of  phase  space  where  the  Kondo  effect  and  RKKY  interactions  do  not  compete 
to  determine  the  nature  of  the  impurity  configuration  at  low  temperatures.  In 
those  regions  of  phase  space  where  these  two  processes  do  compete,  we  expected 
some  novel  behavior  due  to  the  presence  of  a magnetically  frustrated  impurity  con- 
figuration. Specifically,  we  expected  correlations  among  the  conduction  electrons 
induced  by  the  presence  of  a frustrated  impurity-spin  configuration  to  generate 
non- Fermi- liquid  low-temperature  behavior. 

Our  numerical  renormalization-group  calculations  and  the  subsequent  anal- 
ysis have  confirmed  these  expectations  and  validated  our  general  picture  of  the 
interplay  between  the  Kondo  effect  and  the  RKKY  interaction.  For  ferromagnetic 
or  weakly  antiferromagnetic  RKKY  interactions,  our  model  generates  Fermi-liquid 
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spectra.  In  this  region  of  the  phase  diagram  the  low-temperature  properties  are 
consistent  with  a picture  of  three  impurities  being  screened  independently,  irre- 
spective of  any  correlations  which  might  be  induced  by  the  RKKY  interaction. 
Two  exceptions  occur  at  isolated  points  in  the  phase  space  corresponding  to  the 
strongest  possible  ferromagnetic  RKKY  interactions.  At  each  of  these  unstable 
fixed  points  the  three  impurities  are  effectively  locked  into  a spin-|  state.  This 
quartet  is  then  partially  screened  by  one  electron  at  the  R = 0 fixed  point  or  by 
two  electrons  at  the  Jn  fixed  point. 

For  strong  antiferromagnetic  RKKY  coupling  we  observe  two  different  regimes 
of  non- Fermi-liquid  behavior  which  have  no  counterparts  in  the  two-impurity  model. 
One  of  these  regimes  is  destabilized  by  particle-hole  asymmetry  and  hence  is  un- 
likely to  be  reproduced  in  practice.  The  other  (“frustrated-Kondo”)  regime  is 
stable  with  respect  to  particle-hole  asymmetry.  The  low-temperature  spectrum 
and  the  leading  irrelevant  operators  describing  the  approach  to  the  frustrated- 
Kondo  fixed  point  are  unlike  any  seen  previously,  suggesting  the  possibility  of 
novel  non-Fermi-liquid  physics.  We  have  identified  additional  symmetries  present 
at  the  frustrated-Kondo  fixed  point  but  have  been  unable  to  compute  any  proper- 
ties explicitly. 


5.2  Avenues  of  Future  Research 

We  believe  any  future  research  on  our  three- impurity  Kondo  model  should 
focus  on  the  low-temperature  properties  of  the  frustrated-Kondo  non-Fermi-liquid 
fixed  point.  We  feel  our  work  has  identified  the  salient  symmetries  of  the  fixed 
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point  and  that  future  calculations,  perhaps  using  the  techniques  of  conformal  field 
theory,  should  be  able  to  build  on  the  foundation  we  have  established. 

We  believe  that  our  work  also  contributes  to  the  understanding  of  how  ge- 
ometry and,  in  particular,  magnetic  frustration  affects  the  behavior  of  magnetic- 
impurity  systems  at  low  temperatures.  We  hope  that  the  insights  so  gained  may 
eventually  lead  to  progress  on  the  challenging  problem  of  solving  realistic  lattice 
models  of  heavy  fermion  materials. 


APPENDIX  A 

TRANSFORMATION  TO  ONE  DIMENSION 

This  appendix  outlines  the  process  we  use  to  derive  a one-dimensional  rep- 
resentation of  our  three-impurity  Kondo  model.  We  start  with  a more  general 
model  which  describes  a noninteracting  conduction  band  coupled  to  N impurities 
Sj,  where  j = 0, 1, 2, . . . , N - 1 and  N is  arbitrary: 

Hzd  = [ d3k  Y ekcLckCT  + [V'tpl(rj)ipa'(r:j)5(,<(T> 

J O 

-J  Y • &i\-  (A-1) 

We  take  the  conduction  band  to  be  isotropic  in  momentum  space  with  an  energy 
range  ±D.  The  c^’s  annihilate  Bloch  wave  states  which  form  the  conduction  band 
and  represents  the  kinetic  energy  of  these  states.  The  impurity  sites  lie  on  a circle 
of  radius  Rc  in  the  x-y  plane  with  r,  = ( Rc , 7r/2,  <pj)  in  spherical  polar  coordinates, 
The  coupling  parameters  V and  J parameterize  potential  scattering  and  exchange 
interactions  respectively  and  a represents  Pauli  matrices.  The  dimensionless  field 
operators  4>a  couple  the  impurities  and  the  conduction  band,  where 

V><r(r)  = \AV8tt3  J d3kelkrck(T,  (A. 2) 

with  Do  being  the  unit-cell  volume. 
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A.l  N Impurities 

One  identifies  the  conduction  electrons  which  couple  to  the  impurities  by 
first  expanding  the  conduction-electron  operators  and  plane  wave  phase  factor  in 
spherical  harmonics: 


(A.3) 

l,m 

and 

e>k'=4T-£,jl(kr)Yl-J#,mU0,4’)-  (A-4) 

Z,m 

Here  we  use  the  plane  waves  expansion  given  in  Eq.  (16.127)  of  Jackson  [30].  One 
defines  the  momentum  and  position  vectors  k = (fc,0',0')  and  r = ( r,  0,  </ i>)  in 
spherical  coordinates.  The  j/’s  are  spherical  Bessel  functions. 

One  then  uses  the  relation 


Ylm(0,cf>)  = Ylm(01 0)eim*  (A. 5) 

and  the  fact  that  the  impurities  lie  in  the  x-y  plane  {9  — 7r/2)  to  form  the  expansion 
of  the  field  operators, 

Mr j)  = \j2n o/tt  Y,  ilJi(kRc)  £ Yim{ tt/2,  0)e!m^c,w  (A.6) 

l m 

One  substitutes  this  expanded  form  of  the  ipa's  into  the  exchange  interaction  term 
in  Eq.  (A.l), 
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Hint  = -J  {2Qo/tc)  f dk  l dk'^2(-i)lkji(kRc)il  k'ji'(k'Rc)  (A. 7) 

J J w 

x ^y^(7r/2,0)yiw(V2,0)E4w^-'c^'-Eel(m'_m)^  Sr 


Based  on  the  assumption  of  an  isotropic  conduction  band  (ek  — C|k|)>  one 
redefines  the  density  of  states  and  the  electronic  wave  functions  in  terms  of  a 
reduced  energy: 


and 


rlk 

g(s)  = (fto/2ff2)k2- — = q0w2(e) 
ask 


(A.8) 


/ dk 

C-elmo  — W k)  ^ Calmer  i 


(A.9) 


where 


£ = ek/D.  (A. 10) 

Here  the  constant  go  = Q^f)  represents  the  density  of  states  at  the  Fermi  energy. 
One  rewrites  Eq.  (A. 7)  in  terms  of  this  reduced  energy, 


Hini/D  = —AngoJ  I w{e)d£  j w(E')dE'^2(-i)lji(k£Rc)il'ji'(k£>Rc)  (A. 11) 

J J w 

x E n i iw Ee‘{m'~m*‘ si- 


a, a' 


This  can  be  simplified  by  defining 
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S j = Y eij4>/  Sy,  j = 0,  ±1,  ±2,  • • • (A.12) 

y=o 

and 


Cema  — [gm{keRc)\  1 Y<  ii(^'£-Rc)^m(7r/2,  0)  C£/m(T 

l 


(J  — rr)\  I 

= \gm{keRc)]~l  X]iV2/  + lji(keRc)^  ^ + Cdm£r'  (A-13) 


The  normalization  factors  gm{k£Rc)  ensure  that  fermionic  anticommutation  rela- 
tions hold  for  the  new  conduction  electron  operators, 


{c\ma,  C£W}  = S(£  - £')5mm'C',  (A.14) 

which  requires 

°°  ( J 777  ^ I 

M*A)]2  = E(2 1 + 1)U(i.ft)]a|r^|[^r(0)]J-  <A'15) 

In  the  bases  upon  which  the  Sj’s  and  c£ma’s  act,  H^d  reduces  to  the  one 
dimensional  form 


H/D  = [ ds  Y zctmtrCem*  ~ QqJ  [ w(e)de  f w(e')de'  (A.16) 

J m,< j J J 

x Y^,  'Y.  gm{keRc)gm'(ke'Rc)  cl:mcr^o'^(jic£>mic>  • 


This  Hamiltonian  recasts  the  interactions  described  by  the  three-dimensional 
Hamiltonian  given  in  Eq.  (A.l)  in  terms  of  the  interaction  of  an  infinite  number  of 
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channels  of  conduction  electrons,  labeled  by  m,  with  the  impurity  configurations 
acted  on  by  the  Sm. 


A. 2 Three  Impurities 

We  now  show  that  in  the  case  of  three  impurities  placed  at  the  corners  of 
an  equilateral  triangle  in  the  x-y  plane,  the  one-dimensional  form  derived  above 
simplifies  considerably.  In  general  for  a model  describing  N impurities  arranged 
symmetrically  (<£,-  = 2jn/N),  the  Hamiltonian  given  in  Eq.  (A.  16)  reduces  to  one 
describing  N channels  of  conduction  electrons  coupled  to  the  impurities,  where 
the  interaction  term  preserves  the  total  rotational  symmetry  of  the  model,  i.e., 
the  state  of  the  model,  described  by  the  combined  electronic  and  impurity  state, 
retains  the  same  symmetry  under  rotation  before  and  after  all  interactions. 

In  this  case,  we  form  the  impurity  spin  operators  S h from  three  spin-|  oper- 
ators, 


Sh  = £ ei2njh/3  Sj,  h = 0,1,2.  (A.  17) 

j= o 

The  states  acted  on  by  the  S^’s  have  a definite  rotational  symmetry  or  “helicity”  h. 
A helicity-/i  (electronic  or  impurity)  state  acquires  a multiplicative  factor  of  el27r/l/3 
under  a coordinate  rotation  of  27r/3.  Helicity  is  analogous  to  parity  in  the  two- 
impurity  problem  and  reflects  the  three-fold  rotational  symmetry  of  the  impurity 
configuration. 

We  use  the  geometry  of  our  impurity  configuration  to  isolate  the  Y)m’s  which 
contribute  to  the  interaction  term  in  Eq.  (A.  16),  i.e.,  those  which  couple  to  the 
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basis  upon  which  the  S/j’s  act.  Symmetry  permits  just  six  conduction-electron 
wave  functions  of  energy  e to  couple  to  the  impurities.  The  electronic  operators 


Ceha  = 1Xl^/^^(^^/v/3)5^5(m),/iyim(7r/2,0)ceimCT  {h  - 0,1,2) 

l rn 

= [ gh(keR )]_1  ilV2l  + 1 ji(keR/ Vs)  ^ S(m),h> 


\ (f  + raj! 


annihilate  these  wave  functions,  where  (ra)  = ra  mod  3.  Specifically,  c tha  anni- 
hilates an  electronic  state  of  reduced  energy  e,  helicity  h,  and  spin  z-component 
a. 

The  fermionic  anti-commutation  relations  of  the  electronic  operators, 


Ce'h'a' } — E )<W<W, 


(A. 19) 


require  that 


oo  +l  (i mV 

M*)]2  = £(2i  + i)L*(V^)]2  £ iw.'*7rrdr[pr(0)]2. 

1=0  m=—l  ' ' ' 

|[1  + 2sin(x)/:r]  h = 0, 


(A. 20) 


4[1  — sin(a:)/x] 


h = 1,2. 


We  find  the  one-dimensional  form  of  our  three-impurity  Hamiltonian  to  be 


Hw/D  = J deY^Ec\h(Jceha  - g0J  / w{e)de  J w{E')de' 

a,h 

xYlJl9h{k£R)gh'(k£'R)  ^eho\v<rarlCe'h'G'  ' S (h'-h)-  (A-21) 

h,h' 
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Table  1.2  shows  the  basis  we  use  for  the  impurity  configuration.  The  states 
which  form  this  basis  have  a well-defined  total  spin,  2-component  of  spin,  and 
helicity.  In  this  basis  the  one-dimensional  Hamiltonian  conserves  total  helicity 
(modulo  3)  and  is  invariant  under  the  exchange  of  helicity  1 and  2 labels. 


APPENDIX  B 

DETAILS  OF  THE  NUMERICAL  ALGORITHM 


This  appendix  describes  how  one  constructs  a basis  for  and  calculates  the 
matrix  elements  of  the  Hamiltonian  //v+i  as  defined  by  Eq.  (2.29).  Here  the 
procedures  used  are  similar  to  those  used  in  the  one-  and  two-impurity  Kondo 
models  [72,  31,  32], 

One  constructs  //jv+i  from  the  diagonalized  Hamiltonian  H ^ of  the  previous 
iteration  and  the  off-diagonal  matrix  elements  formed  from  /-operators  which  act 
on  site  AM-1,  as  described  in  Eq.  (2.29), 

Hn+ i = YJEm\Em)NN{Em\  (B.l) 

m 

+ £n\P)n+1N+i(P\  + h.C.) |//)  JV+1W+1  (P\ ■ 

h,cr 

One  can  construct  the  basis  for  H^+i  such  that  each  basis  state  |/?) jv+i  has  a 
well-defined  spin,  axial  charge  (isospin),  and  helicity.  Then  the  Hamiltonian  ma- 
trix breaks  down  into  block-diagonal  form,  each  block  containing  matrix  elements 
between  states  with  the  same  total  spin,  axial  charge,  and  helicity. 

When  total  spin  and  axial  charge  are  conserved,  the  process  of  evaluating 
Hjm+ i can  be  greatly  optimized.  In  every  (25  + 1) (2/  + l)-fold  multiplet  of  basis 
states  (or  eigenstates  for  that  matter)  having  total  spin  5 and  axial  charge  / we 
define  the  “top”  state  as  the  one  for  which  Sz  = +5  and  P = —I  [36].  Since  all 
members  of  such  a multiplet  are  degenerate,  we  are  able  to  greatly  reduce  the  size 
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of  the  basis,  and  hence  the  computational  time,  by  retaining  only  the  top  states 
in  a calculation. 

One  forms  each  top  basis  state  \/3)n+i  from  a suitably  symmetrized  product 
of  a top  eigenstate  of  HN  and  an  electronic  state  \1)n+i  composed  solely  of  electrons 
on  site  N + l, 

OO 

\P)n+i  = IOat+i  © \E)n  = s +ns\SE,SE -ns\Sp,Sp) 

ns,nt=0 

X(h,  If  - rii\  Ie,  -Ie  + rii\Ip,  -Ip) 
x|/(ns’ni))jv+1|£(ns’ni)V  (B.2) 

Here 


\l{na'ni))N+i  = N(S+)n°(r)n'\l)N+i  (B.3) 


and 


|E(nj'n<))jv  = N{S~)n‘{I+)ni\E)N,  (B.4) 

where  N ensures  the  normalization  of  each  state.  One  acts  symmetrically  with 
the  spin  and  axial  charge  raising  and  lowering  operators  to  ensure  all  terms  in 
the  expansion  of  each  top  state  have  the  same  spin  and  axial  charge.  The  matrix 
elements  (. . . | . . .)  are  standard  Clebsch-Gordon  coefficients  [58].  In  this  basis,  the 
off-diagonal  matrix  elements  of  H^+i  separate  into  two  components, 


= F(l,l',Sp,Is,P)G(E',E), 


(B.5) 


91 


where  P = N mod  2. 

As  we  shall  see,  F is  determined  once  and  for  all  at  iteration  0.  Only  the 
second  term, 


must  be  evaluated  at  each  iteration.  For  convenience  we  will  henceforth  drop  the 
helicity  label  of  the  electronic  operators  as  the  values  of  the  matrix  elements  do 
not  depend  on  the  helicity  of  the  electrons,  so  long  as  the  quantum  numbers  of  the 
initial  and  final  states  allow  a nonzero  matrix  element. 

Using  Eq.  (B.2),  one  obtains 


G(E,E')  = N{E\f'N(T\E')N, 


(B.6) 


N+liPlfN+lafNvl^N+l  = Sp\Sl,  Sf  ~ 7la;  SE,  SE  + «s) 


x(//3,  -I0\ Ii,  If  - rii ; IE,  -IE  + rii) 


x (//' , If  - n';  Ie1,  —Ie'  + n'| Ip>,  -Ip)  (B.7) 


One  decomposes  the  last  term  in  Eq.  (B.7)  by  noting 


(B.8) 


Thus 
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N 


(£<n-n‘>  1 0 N+i(/(n--)i/i+i./^i^’n;))^i  ® i£'(n>,<)>* 


= (-l) 


(2/5+1) 


N+l 


N+l 


(B.9) 


x *(£("*-n<)|/N,|£'(ns’n<)>N- 


One  evaluates  the  matrix  elements  in  Eq.  (B.9)  with  help  from  the  Wigner- 
Eckart  theorem: 


($/,  Sf,  + n' ; g - n'  + na\Si,  Sf  + na) 
(St',S$;la\Si,S;) 

X (//',  Iii\ 0, 0|/j,  If) 

x n+i(^|/n+i<tIOn+i  (B.10) 


and 


n(E^\  /L  = 

( SE , ■S'jE  ~ n3;  g - n'  + ns\SE',  SE>  - O 
(SE,  Se', ct\Se',  Se>) 

^ ( Ie ) — Ie  + 0>  0|-^£'>  + ni) 

(/£,  — Ie ; 0, 0 1 1 Ef , ~Ie') 

x N(£|/L|£V 


(B.  11) 


Combining  Eqs.  (B.7),  (B.9),  (B.10),  and  (B.ll)  one  can  deduce  F(l,l',Sp,Ip,P) 
by  comparison  with  Eq.  (B.5).  One  simplifies  the  process  of  evaluating  F by 
defining  «(i|/tlO  n in  terms  of  the  iteration-0  electronic  states, 
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»+i(J|/L|Oat+.  = P>l“+W>  offl/llOo,  (B.12) 

where  P = N mod  2 because  of  the  factor  (— 1)N+1  in  the  definition  of  |/)n+i  and 
///  depends  on  the  number  of  / operators  which  make  up  the  state  |/)jv+i-  One 
calculates  F(l,  l',  Sp,  Ip,  P)  once  and  for  all  using  the  iteration-0  electronic  states 
and  the  expansion 


F{1,  l1,  Sp,  Ip,  P)  = 5]  {Si,S[  + ns-,SE,SE-ns\Sp,Sp) 

n3,n's,nj 

x(Sr,  Sf,  + n's\  Se',  Se1  — n's\Sp',  Sp>)  (B.13) 

(Si>,  Sf,  + n's;  \,<J  - n'a  + ns\St,  Sf  + ns) 

x (s,.,s;,-,lc\si,s[) 

(SE,  Se  - ns\\,c  - w'  + ns\SE',  S&  - n'a) 

(Se,  Se',  cr| S&,  Se1) 
x (//,  If  - rii;  Ie,  -Ie  + rii\Ip,  -Ip) 

x(h>,  If,  — nc,  Ie',  —Ie'  + rii\Ip>,  — Ip> ) 

(Ii>,  If  - nj-,0,0\Ii,If  - rij) 

X (It,,  Ifr,  0, 0|/(,  If) 

(Ie,  — Ie  + n,;  0, 0|  Ie1,  —Ie>  + K) 

(Ie,  -Ie',  0, 0| Ie1,  —Ie') 

(J 

One  derives  an  expression  for  G(E,  E')  by  expanding  the  top  eigenstates  of 
iteration  N in  terms  of  the  basis  states  which  form  them, 

\E)n  = ^2cEp\P)n-  (B-14) 

P 
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One  expands  the  matrix  elements  which  form  G(E,  E')  in  terms  of  the  electronic 
states  of  iteration  N and  the  eigenstates  of  iteration  N — 1, 

(£'|/^(T|£',)iV  = CE0cE(3'  N-i(D\  © w(H/jValm,)w  ® 1-^  )w— 1 

w 

= ^2cEpCEp'H(m,m',a,SE',lE',P)SD,D'-  (B.15) 

0,  P 

Here  we  substitute  m for  l and  D for  E to  avoid  confusion.  A relation  similar  to 
Eq.  (B.7)  defines  H, 

OO 

H(m,m' ,cr,  S,  I , P)6d,d'  = (Se,  SslSm,  S^  + ns;  Sp,  Sp  — ns) 

7ls,Tlt 

x(Sm',S^n,  + ns ; So ",  Sp>  - S) 

X (IE,  -IE\Im,  Im  ~ ^ D , + rij) 

X {Jm'  i Im'  - ni'i  Id',  -Id-  + Tli\I , ~I) 
x N{m^’nt)\fUrri,{n,'ni))N  (B.16) 

with  Se(E')  and  Ie(E')  substituted  for  Sp(p>)  and  Ip(p>).  One  transforms  this  into 
the  final  form  with  the  relations 


= 6d,d' 


(B.17) 


and 


S = Se 


S-S*,  = Sjy, 


(B.  18) 
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5 + 5^ 


= 5, 


where  SE  = S'D  as  a result  of  Eq.  (B.17).  Thus 


H(m , ra',  a,  5,  /,  P)  = 

Sn^,ni=o(5  + 5^  — 5^,,  5 + 5^  - 5^/ 1 5m , 5^  + ns;  5 — 5m/,  5 — 5m,  — ns) 
x (Sml,S>m,  + na-,S-Szm„S-  Szm,  - ns|5,  5) 

• (I  + Im  ~~  — 7 T Im1  ~ ^rnl^rni  Im  ~~  n«’  7 — ^m'l  ~7  + ^m1  + W») 


XI 


x (/m/, Izm,  — rii\I  — 7^,,  -7  + 7^,  + n<|/,  -I) 
x 0(m(n*’ni)|/oJrn'(ns’ni))o- 

(B19) 


We  do  not  use  a second  set  of  raising  and  lowering  indices  (n's  or  n')  in  this  case, 
because  of  the  equivalence  of  the  eigenstates  of  iteration  N - 1.  One  finds, 


G(E,  E')  = J2  cEpcEp'  77(m, m' > a>  Ie '*  p)^d,d',  (B-20) 

P a 

which  requires  knowledge  of  the  coefficients  cEp  which  go  into  the  expansion  of 
\E)n  at  each  iteration  as  well  as  a knowledge  of  how  the  basis  states  \(3)E  were 


constructed. 
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